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TECHNICAL P^EMORANDUM 


EQUILIBRIUM AND STABILITY OF A SATELLITE INFLUENCED 
BY GRAVITATIONAL AND AERODYNAMIC TORQUES 

INTRODUCTION 


The stability of satellites under the influence of only gravity gra- 
dient torques has been treated by several investigators [1-4]. The inclu - 
sion of the effects of aerodynamic torques was considered by Nurre [5] 
and Sperling [6], In addition, the instability produced when a gravity- 
gradient stabilized body is subjected to a torque which forces it slightly 
away from the true gravity- gradient equilibrium position (i.e., axis of 
minimum moment of inertia aligned with the radius vector and the axis of 
maximum moment of inertia aligned with the perpendicular to the orbit 
plane) has been studied by Garber [7]. The purpose of this report is 
to document a method which can be used to determine the equilibrium con- 
ditions and investigate the stability of these conditions for satellites in 
circular orbits under the influence of gravity- gradient and aerodynamic 
torques. The procedures considered are demonstrated using the former 
Sky lab satellite as an example. 

The satellite is assumed to be rigid and the atmosphere constant for 
a given altitude. A linearized infinitesimal analysis is used to assess the 
stability of cquilibi’ium positions found by a variation of the Newton- 
Raphson method. 


COORDINATE SYSTEMS 
Orbital Reference System 

The orbital reference system has its origin at the center of mass of 
the satellite and I’otates uniformly at the rate of one revolution per orbit 
(n). The axis is positive in the direction of flight and lies in the 

orbit plane. The axis coincides with the orbital radius vector and is 

positive away from the Earth. The axis is parallel to the orbit normal 

and the orbital angular momentum vector. 


Body-Fixed Systems 


Two body fixed systems are portrayed in Figure 1. The Xg, Yg, 

Zg is a geometric system located at the center of mass. The principal 

axes system, Xp, Yp, Zp, also has its origin at the center of mass and 

is oriented such that the product of inertia terms vanish. It is in the 
principal axes system that the equation of motion is written. The rota- 
tion from the orbital system to the principal axes system is accomplished 
by three positive rotations ij/g, and about the three successive 

axes 1, 2, 3, respectively. Thus, 

P i. ^ 0 


where over-scored symbols represent vectors 
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- Si Sj Sj 

S, 0, + Cj 




In these matrices and s^ represent the cosine and sine of tpj^, respec- 
tively, etc. 



EQUATION OF MOTION 


The equation of rotational motion of the satellite in the prineipal 
axes system is 


Ivv* + w'* X Ivv* “ E 


where 


I “ the constant moment of inertia matrix 


w* s angular velocity of the satellite with respect to inertial space, 
represented in the principal axes system 


L = external torques, wluch for our purposes are gravitational and 
aerodynamic toi’ques, 


In particular, wo have 

w * = w + n , 


where 


3 


w != angular velocity of the body system with respect to the orbital 
system 


Q s [i|;] is the angular velocity of the orbital system with 

respect to inertial space. For representation in the body sys- 
tem it is multiplied by the transformation matrix 


Now, 
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Therefore , 

w * = w + ?l s w + ttl + IiIj] 

In our case n’t* = 0 and 
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Thus the equation of motion is 


I (vv - w X Sl) + (w + u) X I (w + II) - X = If . (1) 

« 

EXTERNAL TORQUES 


The gravity gradient torque on a rigid body in circular orbit about 
its three principal axes is: 


= 3 ir r X U. , 


Where 


r X ir = 


(^3- V ^^1^2 «3> 

(Ij^ - Ig) Cj^ Cg (Sj^ Sg “ Cj^ Sg Cg) 


(I2 “ Ip Sg - Cj Sg Cg) (Sj^Cg + Cj^ Sg Sg)^ 

m 

The aerodynamic torque is not so easily written. It is of the form 


L = K 
a 


L, 


2 

where K = l/2pADV in which p is the atmospheric density, A is the 
reference area, D is the reference length, and V is the magnitude of the 
velocity. 


For i = 1, 2, 3 
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6 

E Oa) 

C®1 

6 

+ E ’*' ®®ijk ®*” 0“) 

k""X 

in wliich AA, AB, BA, and BB are the coefficients of the double Fourier 
series representing the aerodynamic torque in the principal axes system. 
The angles a and ^ are the total angle of attack and the aerodynamic roll 
angle, respectively. These angles are' depicted in Figure 2 and are 
dollned as follows: 




= cos”^ (Vj^) 
i|) s tan“^ ’ 

where Vj^, V 2 , Vg are unit vectors of the velocity vector expressed in the 

geometric body reference frame. It is in this reference frame that the 
aerodynamic moment and force coefficients are rendered [ 8] . 

EQUILIBRIUM 

The equilibria of the linearized equation of motion (1) are solutions 
of 


V = i' 2 X I U - L = 0 , 


To find an equilibrium point it is necessary to find a zero of v. However , 
the aei'odynamic torque component of the external torque is si^ficiently 
complicated to preclude the analytic determination of zeros of v. Now, 

V is a function of 



6 




and on iterative procedure like that of the Newtou-Rophson method can 
be used to solve for the values f which render v » 0. We then have 



where succession approximations of v are noted by subscripts i and i+1. 
Then for to be zero, it is necessary that 



LINEAR STABILITY 


For a linear 'stability analysis, wo let i|j = tji* + ^ where Iji* 'is a 
zero of V and 



is a small displacement from an oquilibriuoi point. Note that ^ ~ ^ and 

• • • t 

ip s ■ 3 . The equation of motion [equation (1)] can be expanded to 
Iw “ I(w X W) 4- w X iw + li X iw + w X ivv + V = 

• • 

After linearizing and x’ocalling that w was of the form [Q] ip and w of 

• • 

of the form [Q] ip, it is possible to arrive at 
[A] 6 + [BJ ^ + [C] « = 0 
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or, 


t + |AJ“^ IBI 5^ + [CJ ?T a H 


where 


lAJ s tl] tQ] , 


[B] = 


(I3 “ Ij “ 1 2^ ^^3 


(I3 - 13) ^3 -dg - I3 - Hg 

-(I, - lo “ I«) «, 


1 

0 


(1 2 - I3 - I|) Rg “C*! “ *2 " * 3 ^ % 


1 ‘2 “ * 3 ^ 'd 

0 


p - iL}? 


where I. and are the components of I and R, respectively. (The matrix 
A can become singular when the cosine of i/zg equals 0 .) If we let 


then 


It is then possible to write equation ( 1 ) in terms of c as 


c + 


^ -U 


a‘*^c a“^b 


c= 0 


( 2 ) 




where 


and 


[ 0 0 Ol 
0 0 0 
0 0 oj 


U = 


10 0 
0 10 
0 0 1 


To determine how unatable the equilibrium point is, it is necessary to 
obtJiin the oigenvoluos of equation (2) and compute the unstable time con- 
stant or the logarithm of the amplitude over one orbit (logarithmic 
increment) . 

A computer program was doyelopcd to first determine coaibinations 
of \j> 2 , '/'g which would render v - 0 and then to determine the eigen- 
values of the system for each such attitude. In the case of the Skylab, 
it can be demonstrated that v is not a unimodal function so that finding 
all equilibrium attitudes cannot be easily done except perhaps by exhaus- 
tive search. However, for purposes of demonstrating the operation of 
tliis program, an initial starting sot of angles was chosen close to an 
equilibrium point and after a fixed number of iterations, the stability of 
the Skylab for that equilibiauni point was assessed for combinations of 
altitude and density. 


COMPUTER PROGRAM 


The form of the program is relatively straight forward. Input 
quantities are read in first. The aerodynamic quantities, which are the 
Fourier coefficients for aerodynamic forces and moments in the geometric, 
system, are transformed into the Fourier coefficients for ilm aerodynamic 
moment in the pidncipal axes system. This is done in Subi’outine ATRAN. 
The progranj then accepts combinations of density and altitude and initial 
angles i{), <*, ij)^^ for starting the iterative search for combinations of 

iljgi ipg (Euler angles) such that v = 0, Upon reaching the iteration limit, 

the program outputs the values of ij/g, ti/g in the form of t|), a, 4>g 

(for ease of visualization). The magnitude and components of v are also 
displayed for visual determination of how close to equilibrium the last 
angles are. These quantities, as well as the associated eigenvalues, are 
then stored for possible future plotting. 


XO 


A listing of the main program and associated subroutines is given 
in Appendix A. A listing of major program variables and input data are 
supplied in Appendices B and C, respectively. In Appendix I) are 
presented the results of a sample cose and some relevant plots. 

The program was written for the Honeywell CP-V System. It should 
be noted that there are two machine dependent parameters; one in Sub- 
routine BALANC and the other in Subroutine HQR2, which must be altered 
for use of this program on other computing systems. 
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PB0CR4M TO FIND POUILIBRTUM ATTITUOfS FO« * 8*TEUITK ()NPF» 
CUAVITT GRAOlfNT AND ir»OOVN4Mir TO»QUF8 

FOURIER COFFFlCirNTS FOR 3EROPVNRM1C EFFECTS FRF INPUT IN THE 
FORM GIVEN RV NORTHROP “EMO m- 9?50«75-I R7 BY 
P.F, r.VORFI 


DIMENSION Tr.P(5,S).TPr.(3,3KPM.‘S(3,3),nNIJS(3#3),PNII8I(3.5). 
lf)ELP8I(3),PA8f3),PSH3).Von',f' !I0(3>,C(3)»S(3),V(3),*NU(5)» 
2FN(S)»PAMA(3),p»mp(S),pp 8(7),^ AtiS(5,^),XI(3),PFSH3»3) 
dimension ALAn( 6 . 6 ),*LR 0 ( 6 .t >, '*0f6)#*l»(6.6f6).AL8(6.6»6I. 
1ftLF(6,6,6),ftLH(6,6,6),PSlG(i),win),MH3),PPHI(3),APHIf5,3) 
DIMENSION AHAn(3,2).»K80(3,?>,*KA(3,?,6),*KB(3,?,6)# 
lRK4(3,?,6),»KH(3,2,6),4DtiM(3,31,TMM(6,6).r)ELPG(3),TC(6) 
OIMFNSfOM T,R(6).6t (M,Z(6,6),FV(6) ,IV(6),OltP(3,3),EUM(3.3) 
DIMENSION OI(3.3>.B(3,3),Pr.r.S(3.3),F»(6),7ET*(6),llf3,6) 

REAL PI/3. l4tS9?653590/ 

PAO=S7. 298779515 

INPUTS 

PrFFRfl.Cf DIA^FTfR IN FFFT 
PFrFPFNri 4PFA IW FFETP 
» U^*«{ H PF ITFRATin-S 
. TOP 

►>'nMFUTS OF inertia AftniiT PRINCIPAL AYES IN SLUR-FEET2 
ifiCATinN nr CG n gfomfipic rnnwoiNATES in feet 

HEAn(5,700)nPEF,APEF,LIMII 

Rr 4 i){ 5 , 70 | )f fTGPFL,R),M*) ,TI,L *1 » 3 ),yn ,*I 2 ,yi 3 ,)(CG,YCP,,ZCR 
PEAnf 5 , 702 )Mi ,N?,M 3 
Rni)f 5 , 703 )(Aftn(i),ls| .Mil 

no 5 1 * 1 . n 

RtAn(5.70 3) (Al AP(T,,n,Js1 ,^V) 
s RFAn(5,703) (ALftOd . Jl.Jsl.M?) 

no 7 i«i,vi 
no 6 jsi ,M? 

RFAniS, 70 3 )(ALA(I. J.K),K=I .'-T) 

6 REAP(5,7U3) (HLA(I,J,K),K=1 ,"3) 
no 7 .1=1,"? 

REAOf5,7niHALH(I,J.*<) ,K*1 ,*'31 

7 REAn(5,703Uftlft(I(.J,t<),K = l,M8) 

700 FORMAT(?r,/n 

701 F0R'’AT(3G) 

702 FPR"ATf3T51 
705 FORMATFaF |?.0) 

RF»?09?57?1 .7ft5 
GNsl .4fi765 39jri6 
on ft L*l ,3 
no ft M»i ,5 

ft TPR<L,M)sTRPfNj ) 
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c 

C INTEOIA TERMS 
C 

XKDOMtl 

KI(2)«*12 

XI(3)««I3 

XMUX13-X12 

)(M2»)m-XI5 

XM5xWl?-)(ll 

KK1*XM1/XU 

X«?«XM2/XI2 

*rii''cr./oPKF 

VD«VfC/r>REF 

zd*zcg/ohff 

c 

C CONVFROIflf; TEPMS FOR CHANGTHR AEROnvNAM|C TOROtlF 


TO 

PRTNCIPAI 

AXF 

S 

reference 

FRAMF 

no 

0 1=1.3 





■v)(l 

.1 )»Z0MPG 

(I. 

?1 

-Yn*TPr,(T, 

>i) 

U(I 

.2Y»Yri*TPG 

(I. 

1) 

•XD*TPR(I, 

»2) 

un 

,3)*ZD*TPG 

(1. 

11 

-X0*TPG(1, 

>31 

i|(T 

.0)*TPG(T, 

2) 




IHl 

,5)*IPG(I. 

31 




11(1 

,6)«TPG(I, 

11 





0 rONTlNUt 


C CONVERSION. OF AEROOYNAMIC COEFFICIENTS FROM r.fnMFTRlc 
C RFFERFMCF FRSMf to PRINCIPAL AXFS REFFRNCE FPAMf 
C 

CALL ATkANtALAO, ALR0»AI A, ALH.61 i i PLR. AK AO, AKHO, AK A, AKR, 
1RK»,BKP.U1 
C 

r PHOF « ATN'OSPHFRIC DENSITY IN KG/M-J 
C (PROGRAM TEPMifjATED IF RHOK t 0) 

C ALT ■ ALTITIJDE in NAUTICAI MUfS 
C 

lOnn REAO(5.70t 1RH0K,ALT 

TF(RHOKl1001»100?,in01 
1001 on 10 1=1.3 
no 10 j«i.6 
10 TMMn,J)«0. 

TMMM ,u)s-1 . 

TMM(3,6)«-1 . 

RH0*RHrK/S.1537OFF2 
ROPP*RF*ALT*6.0761 1506E+3 
nMC»SOPT(r.‘VRnRB)/RORB 
nM?sOMR«OMr, 

GG*5.*nMR*nMG 

OENeRHO*ORFF*ARFF 


IQS. 

106. 

107. 

10 «. 

109. 

no. 

in. 

112 . 

m. 

no. 

ns. 

116. 

117. 

no. 

119. 

120 . 
121 . 
122. 
123, 
120 . 

125. 

126 . 
127. 
120 . 

129. 

130. 

131 . 

132. 

133. 
130, 

135. 

136. 

137. 
136. 
139. 
100 , 
101 . 
102 . 
103. 
100 . 

105. 

106. 
107. 
ion. 

109. 
ISO. 
151 . 

152. 

153. 
150. 

155. 

156. 


ClOa.SnOENnVMAGoVMAG 

9EAD(Si7OnO870 

C 

c initial ANGIESi 

f 

C PSIOCn • ALPHA « TOTAL ANGLE OF ATTACK 
C P810(2) * PHI ■ «CU ANCLF 

C PSI0I3) • PHI HETA I BONK ANGLE 

C 

C CONVEPSION TO EULER ANGLES PSI(1)» P9I(2), PSH3) 

DO n LJ«1»3 

n prhi(lj)«psio(lji/rao 

cpi«cos(ppMi(m 

CP2«Cns(PPHi(2)) 

CPl»C0S(PPHT{3n 

sPi»8rN(ppMinn 

SP2«SIN(«»PN|f2)) 

SP3»S1MPPHI(3)) 

APHI(l,n« rP2 
APHI(n2)i-SP2*SP3 

APHI(|,3)i SP2*CP3 

APHIf2,|)» SP2‘SP1 

APHI(2.?)*.CPl*CP37SPnrP2nSP3 

APHH?,3)»-rPi*sP3-SPnr.P2*fP3 

APHl(3,n» SP2*CP1 

APH!C3,2)« SP1*CP37CP1*CP2*SP5 

APH 1 (H, 3 )« SPnSP5-CPl*CP2Ar.P3 

CALI. MAT3(TPfwAPHI,ARllM) 

PST(n«ATAN?(-A0UM(3»21»Api|M(3,5)) 

PS! (3)»ATAN2(-ADIIM(2, n.AmiNd , j)l 
ADitnA0l(V(3,5)/r0S(PSt(l )) 

PSI t2)«ATAN2(AmiM(5.1 ),AfHl) 

C 

r BEGIN ttepation loop 

C 

nn too KotiNT«l .limit 
ci«cns(P5i(in 
C2«cn8(psi(2n 
C3sCnS(PST(3)) 

Sl«v5IN(PSI(l)) 

S2«SIN(PSI(2)) 

S3*SINfPSI(3n 

csi«n*ri-si*st 

CS2=C2*r?-S2*S2 
CS3»C3*C3-S3«S3 
SC3I*3.*C1 *C1-S1*S1 
■ni(l,n»C3/C2 
0I(l,2)s-S3/f? 

01(1,3)30. 

OK?,! )tS3 
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185. 
180. 

185. 

186. 

187. 

188. 
1"9. 

190. 

191. 

192. 
195. 
19«. 

195. 

196. 

197. 

198. 

199. 

200. 
201 . 
202. 
205 . 
200. 

205. 

206. 

207. 

208. 


R1(2,2)»C3 

01(2*3)a0. 

01(5,1)«-S2*C5/C2 

am»?)»82*S3/C2 

0I(3»S)«1. 

B(l,l)»OMC*(OfKl8C82-l.)*Sl*85t(XKHl,)*Cl*82*C5) 

B(1,2)«-0MG*(XK|-1,)*S|*C28C5 

8(l,3)«nMr,i(XKl*1,5*(Cl*C5«Sl *82883) 

«(2.n«-0‘’G*((KK28CS281,)*Sl8Cl-(XK2-l.)8Cl*8?8S5) 

9(2,2)«-OMG.(XK?*l.)*Sl8C2883 

8(?«3)80MC*(XK2.1.)*(C18SS88I8S28CS) 

11(3,1 )«0MG*((XK5*CSS-1.)*C1*C2-2.*XK58S18828C28838C3) 
8 ( 5 , 2 )inMr.*((XK 58 CS 58l.)*Sl*8282.8XK38Cl8S38C5) 

P(3,5)«0. 

C 

C PfiGS IS THF PARTIAL Of OMJgA K I8nMfGA - GRAVITY GRADIENT 

r WITH RrSPFCT TO PST AND FXPK'FSSfD IN TNf PRINCIPAL AXES 

C RfFERENTF FRANF 

r. 

PGGS(l,l)»-0.*n^2*C^•YMl.(C3•CSl-2.*Sl8Cl8S^•S3) 
PGr.fi(l,21s-)fW) *oV28(CS?*S3*SC31-<I.*81*C1*S2*C5) 

PCGSd ,3)»-0M2*C2*XM18(SC31*S28C3-0,*S1*C18S3) 
pr,GS(2,n«-o.*OM2*c?*xM?.{838Csi*2.*si8ri*s?*r3) 

Pr.GS(2,2)a XN?*0^*?8(CS2#C5*SC518«.8S1*C18S2*83) 

PGGS (2,5) 5-nM2*C2*XM?. (SCSI ♦S28S5*0.*S1*C1*C3) 

. POr.SfS.ns u.*nMp.XM3*(CSl8S2*CS3-2.*Rt*Cl*(1.8S28S2)8S58r.5) 
PGCS(3.2)* 2.*0M2*C28XN3*(SC3l8v82«S38C3T2.8.Sl*Ci8CS3) 
PGGS(5,5)* nM?.XM3*(rS3*(S?882*SC3l*Cl*n-3.*31*SI )- 
1168Sl*Cl8S?*S3*r5) 

C 

C UNIT VFLOriTY mFCTOR IN PRINCIPAL AxfR RFFERF^CF FRA«E 
C 

V0(n«C28r3 

V0(2)*-C2*S5 

V0n)»S2 

C 

C TRAMSFt'R^^ATiriN OF UNIT YFLOflTV VECTOR TNTO GFONFTRIC 

C REFERFNCF FRAVF 

C 

CALL M«TVEr(TnP,VO,V) 

VM )*AMAXf A'*J,’( 1 ,o,V(l) ),-1 ,0) 

c 

C ALPHA S TOTAL ANGLE OF ATTACK 

C PHI s ROLL ANGIF 

C 
C 

Al.PHAsACns(Vf 1 ) ) 

PHI=ATAN2(V(2),V(3)) 

CONSTlsl .n-V(l )*V(1 ) 

IF(CnNSTl)20,25,?0 
20 rONRT2»l ./SORTCCONSTl) 

GO TO 30 
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20«, 

2 t 0 . 

2tt. 

212 . 

2IS. 

2U. 

21 5. 

216. 
217. 
21 «. 
21 «. 
220 . 
22 1 . 
222 . 
?? 1 . 
22 «. 
22^. 
226. 
227. 
? 2 «. 
220 . 
250. 
231. 
252. 
255. 

254. 

255. 

256. 

257. 
256. 
250. 
240. 

?ai. 

242. 

24 5. 

244. 

245. 

246. 

247. 

248. 

249. 

250. 

251. 

252. 
255. 

254. 

255. 

256. 

257. 

258. 
250. 
260. 


25 C0HST2«l.5f70 
50 CONTlNUr 

c 

C C8LCtlL*TinN or A|(»00V»t4MIC TOUmif 660 P*»T!AtS OF 
C 6CR00V»‘*M1C TOROlir WITH RfSPFCT TO 6LFH* «N0 PHI 
C 

CALL 26CPO(«)(AO,«KIIO,6)(*««FR,RK«,nRR.«lPH«,PHT,CLA« 

IPUHA.PAMPpl^'l 

C 

C PAS * PARTIAL OF ALPHA pTTH RKSPtCT TO P8I 

C PPS • PARTIAL OF PHI PITH RfSPICT TO PSI 

C 

PA 8 (n« 0 . 

PAS(2)»f S2*f TPCn «1 )*C5-TPr.(2,U*SS)-TPr.(5»n*C2I*C0H8T2 

PA8(5)»C?«(TPG(l,n*8|*TPC(2.U«C5>*C0NST? 

ppsniao, 

PPS(2>«<V(2)*(S?*(TPC(|,5)*C5-TPf.(2,S)*S5l-TPr.(5,8MC?)- 

IV{5)*(S2*(TPf,(|,2)»f.5.TPr.(?,?)*85)-TPC(3.2)*C?>)/COHSTl 

PPS(3)«C2MV(?)*(TPCn,5)*S5tTPC{2,5)*C5)-V(5)*(TPC 

m, 2)*85*TPr.(2.2)*C5))/fOK'STl 
C 

C CALCULATIC'N OF PARTIAL Of AFROPVNAHIC 

C TOROIIE HlTH PtSPECT TC P81 

r 

no 40 1«1.5 

00 44'5«1,5 

40 PAMS(l,J)»PAMA(I>*PAS(.|)4P»MPn)*PPS(J) 

c 

C CALCULATION OF PARTIAL OF FXTFRNAL TORQUF 

C WITH RFSPFCT TO PSI 

C 

on 50 I»l.3 
DO 50 J«1,A 

pNusd, n«p(!osn, j)-p*Hs(i»j) 
pA8tn,j)»PMisn,j)/Ain) 

50 (3NllS(dJ)«PNUSn.J) 

(; 

c calculation of nu 

c 

)(NtHn»-nM2i*Mt *C2* (4.*si *C1 •C5 aS2*S3*SC51 )-AMU ) 

XNU(21e 0A'2*»tH?»r2*(S?*C5*SC5l-4,*Sl*CI«851-AM(2> 
XNU(5)sO‘A?*XH5*(U,tSl«Cl*S2*CS5A85*C5*(52*S2*SC5UCI*Ct- 

n. *Sl»vSn)-Ah(5) 

CALL minV(0NIIS.1»0I.11,MT) 

C 

C PNUSI ■ INVFRSF OF PARTIAL OF NU MTH RESPFCT TO P5t 

C 

00 60 1*1,5 
DO 60 J«1 ,5 

60 PNUSI(l,J)»OMt'S(I,J) 

CALI MAtVEC(PNllSI,XNU,nFLPSn 
XNI|MAG«0, 
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Hi, 



00 70 I«l,5 

Hi, 



XNU^Ae«XNUMAO*XNU(p*KNUn) 

?6S. 



DEU»G(n«-nCLPSHI)*960 

Htt, 



PSIC(I)iPSHn*P«D 

Hi, 


70 

0ELP8I(!)»-0ILPSim 

Hb, 



KNlIMftCvSMRf (XNUPAG) 

H7, 



ALPHAG9ALPH«*HA0 

Hi, 



phic*pn:*pad 

Hb, 



MXM«.(TCP(l,l)*(CMSSA91*82*C3)tTGP(l»2)*(Cl*C!-9l*S2*85)- 

?70, 



1TCP(1»3)*8I*C2) 

271. 



UXO«Tr,P(I,n*(Sl983-CI«82*C3)tTGP(l,2)* 

272. 



IC8l9C37n*52«S3)*TCP(l»5»9Cl*C2 

273. 



PHrBC«ATAN2nJX*l,UXO)*RAn 

27«. 



90 90 !.«|,3 

275. 



P8HU9AMOO(P8I(U7DELP8I(L),2,0*PI) 

iH, 


90 

lF(A98(P8l(U).fiT.PI4.0nPSI(l )«P8I(L)-8IGN(2,0*PT»P8ia) J 

ill. 


9ft0 

F0RMATf//,5X»9HALmtl0f ,P6. 1 , 3X.9HOE-MSITY ,E11.5) 

27«. 


901 

F0R«AT(?X,<IHPHI ,FI1,6,2X,6MALRHA ,FI0,6,?X,5HPH!8 »F10.6) 

279. 


902 

FORMAT(lX,SF15.0»IX.El5,9,y) 

2S0. 


903 

format ( 7X, «H»r AL, l5K,9H|MAf.,13X#9HFRFO#5X.9H2ETA,6X,2hrC) 

291. 


909 

FO»MATOX,M6,a,lX,fl6.0,2X,FlO.5,2K,r6.O,2X#E9.2) 

29?. 


100 

CONTINUE 

295. 



CALL W4T5(0T»PASI,0U^') 

PH, 



CALL MAT5(01,n,ElJM) 

295, 



00 110 I«1,3 

206. 



J»T*S 

207. 



no no (‘■ 1,3 

299. 



LX«X*3 

209. 



TMM(J,»<)BnUM(j,K) 

290. 


110 

TMM{J,LK)«El'**(I»t') 

291. 

r 



292, 

c 


nETF9Mf'''ATlON PF FICrNYALUFS 

295. 

c 



29U, 



TFRR*0 

295. 



CALL EIS(TMM,*^■R,t• l,Z,FVnv»6,6,IFRR) 

296. 

c 



297. 

c 


Of TF.RMj\ATinM Of unstable TIME CnUSTANTS 

299. 

c 


OR LOGARITHMIC INCREMENTS 

299. 

c 



300. 



CALL E!SnRT(^R,Wl,FR,TC,7FTA,0Mn,Pl,6) 

501 . 

c 



302. 

c 


OUTPUT TO TFRminAL SCREEN 

305 . 

c 



3oa. 



*<(RITEnC2,90n)ALT,RHOK 

305. 



rtRITE(lf'2»90|)RHTn,Al PHAG,PMIHG 

306. 



WRITEn02,9n2)yNU,XNUF4G 

307. 



rtRITE(|n2,903) 

309. 



wHITE( 10?,909HWR(T),1»I(I) ,FP(n,7ETAm,TC(n,T*1 .61 

309. 

, c 



310. 

c 



311. 

c 


•WRITE miTruT filf for plotting 

312. 

c 



313. 



wRITEn01RMnK,ALT,WP(n,WI(l ).WR(3),WI(3),WR(5),WI{5),PHTG 

31". 



1 ALPHAG.PHTRG.XNiJHAG.TCd ).TC(3) .TC (51 

315. 

c 



316. 

c 


PAUSE for terminal SCREEN RESET 

317. 

c 



319. 

c 


PAUSE 

319. 

c 



120. 



GO TO 1000 

321. 

10^2 

STOP 

322. 



fnd 
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1. 



SUBROUTINE MINV(a,N,D,L***I 




MIN 

Root 

2. 

C« 






*MIN 

0002 

3. 

C 






MIN 

OOOS 

«. 

c 


MATRIX INVERSE* OF AN NXN MATRIX-A INTO 

THE SAME MATRIX-A 

MIN 

OOOR 

5 . 

c 






MIN 

OOOS 

6. 

c 


-1 0 m OUTRUT 

determinant Of ORIGINAL A 

MIN 

OOOR 

7 . 

c 


A ■ A L* M ARE 

THD HORXINC INTEGER VECTORS 

MIN 

0007 

a. 

c 


OIMENSIONfO 

IN MAIN* OE LENGTH 1 

N*»IN 

OOOB 

9 . 

c 






MIN 

0009 

to. 

c 


NOTE THAT THE ORTfilNAL A-maTRIX TS DESTROYED 


MIN 

0010 

n. 

c 


ALSO MOTE THAT A CHEER ON D SHOULD BE 

MADE IN 

CASE MATRIX-A ever 

MIN 

OOtl 

12. 

c 





BECOMES SINGULA* 

MIN 

0012 

IS. 

c 






MIN 

0013 

14. 

c* 



»*« 



AMIN 

0014 

15. 



DIHENSfON A(l )*(.(n*M(|) 




MIN 

0015 

16. 

c 


SEARCH FOR IAHGEST ELFHENT 




MIN 

0016 

17 . 



.0 




MIN 

0017 

la. 



NK«-N 




MIN 

0010 

19. 



DO ao K«i *u 




MIN 

0019 

20, 







MIN 

0020 

21. 



L(K>»< 




MIN 

0021 

22. 



M(X)IK 




MIN 

0022 

25, 



KK»MK4K 




MIN 

0023 

?«.’ 



RIGAbA(KK) 




MIN 

0024 

25 . 



no 20 jsK,f) 




MIN 

0025 

26. 



iz*N*(j-n 




MIN 

0026 

27. 



on 20 !«K,fj 




MIN 

0027 

2 a. 



1J«IZ*T 




MIN 

0028 

29. 



IF( ARS(BIRA)» AflS(A(IJ))) 15,20.20 




MIN 

0029 

30. 


15 

BIGAaA(IJ) 




MIN 

0050 

51 . 



L(K)*I 




MIN 

0031 

52. 



MUlaJ 




MIN 

0032 

53. 


2ft 

CONTINUE 




MIN 

0055 

34. 

c 


INTFRCHANOF ROftS 




MIN 

0034 

35. 



JsL(K) 




MIN 

0035 

36. 



IF(J-H) 35,55,25 




MIN 

0036 

37 . 


25 

KlaX-N 




MIN 

0037 

3 a. 



DO 50 Tsl,N 




MIN 

0038 

39. 



KI«K I4N 




MIN 

0059 

40 . 



H 0 L 0 a-«(an 




MIN 

0040 

41. 



JlaRI-R+J 




MIN 

0041 

42. 



A(KI)sAfjn 




MIN 

0042 

43. 


50 

A(jn *HniD 




MIN 

0043 

44, 

r 


TNTERCHANOF columns 




MIN 

0044 

45 . 


55 





MIN 

0045 

46, 



IF(I-K) 45,45,38 




MIN 

0046 

47 . 



JP«N*(T-t ) 




MIN 

0047 

4 fl. 



no 40 j»i,r' 




MIN 

0048 

4«. 



JK«NK+J ' 




MIN 

0049 

50 . 



JJsJPf J 




MIN 

0050 

51 . 



HnLOa«A(.rK) 




MIN 

0051 

52. 



A(JK)«A(JI) 




MIN 

0052 


20 


53, 


60 

AtJI) tHOLO 

MIN 0053 

54. 

c 


OTVIPF CflLUMN BY MINUS 8IV07 (VALUE nr FlVn? CLCWfNT 18 

MIN 0056 

55. 

c 


COMTAIfEfl IN HIGA) 

MIN 0055 

56. 


65 

ir(AHS(8!CA).Uf-20)66»66.68 

MIN 0056 

57, 


66 

DIO.O 

MIN 0057 

58. 



RETURN 

MIN 0058 

56, 


68 

DO 55 I«1,N 

MIN 0059 

60. 



ir(I-Kl 50.55#50 

MIN n060 

61. 


50 

IKbNK^I 

MIN 0061 

62. 



A(IK)*A(IiO/(-Pir.A) 

MIN 0062 

63. 


55 

CONTINUk 

MIN 0063 

66. 

c 


REDUCE M8TRIX 

MIN 0066 

65. 



DO 65 I»1,N 

MIN 0065 

66. 



IK«NK»! 

MIN 0066 

67. 



HOLDiAdlO 

MIN 0067 

68. 




MIN 0068 

60. 



PO 65 J«1,N 

MIN 0069 

70. 




MIN 0070 

7t. 



1F(1-K1 60#65»60 

MIN 0071 

72. 


6ft 

1F(J-V1 6?»65,62 

MIN 0072 

75. 


62 

8J9IJ-1+8 

MIN 0071 

7tt. 



*nj)«HPLr>«ft(K.i)»A(lJ) 

MIN 0076 

75. 


65 

conttnuf 

MIM 0075 

76 . 

r 


OIVIDF ROW nv PIVOT 

MIN 0076 

77. 



K.J«8-N 

MIN 0077 

78. 



DO 75 .1*1. M 

MIN 0078 

79. 



KJ«KJAM 

MIN 0079 

«0. 



IF(I-K) 70,75.70 

MIN 0080 

81. 


70 

A(XJl«*(K,n/HIGA 

MIN ft08l 

82. 


75 

CONTINUfc 

MIN 0082 

83. 

c 


PROmiCT OF PIVnTS 

MIN 0083 

86. 



0*D*H!r,A 

MIN 0086 

«5. 

c 


REPlACe PIVPT BY RFCTPROCAL 

MIN 0085 

86. 



A(KK)«1 .0/PlGA 

MIN 0086 

87. 


80 

CUMTINWE 

MIN 0087 

88. 

c 


final R'jvh AMD column IMTERCHANGF 

MIN 0088 

89. 



K*N 

MIN 0089 

90. 


100 

K*(K-n 

MIN 0090 

91 . 



IF(K) I50.1'50,105 

MIN 0091 

92. 


105 

I*L(K) 

MIN 0092 

93. 



1F(1-K1 120,120,108 

MlN 0093 

96. 


108 

jn«N«(K-l ) 

MIN 0096 

95. 



rH*N«n-i ) 

MIN 0095 

96. 



DO no .1*1 ,N 

MIN 0096 

97. 




MIN 0097 

9fl. 



HnLP*Af.lh1 

MIN 0098 

99. 



Jl*JR*J 

MIN 0099 

ion. 



A(JK1«-A(JT) 

MIN 0100 

101. 


1 10 

A-(jn sHDUD 

MIN 0101 

102. 


120 

JsM(K) 

MIN 0102 

105. 



IF(J-K) 100,100,125 

MIN 0103 

106. 


125 

KT*K-N 

MIN 0106 

1 05. 



DO 130 !»I,N 

MIN 0105 

106. 



KI*8UM 

MIN 0106 

1 07. 



HOLn*ACMl 

MIN 0107 

108. 




MIN 0108 

109. 



A(Kn*-A{JI) 

MIN 0109 

no. 


130 

A(JI) *HnLD 

MIN 0110 

111. 



GO tn joo 

MIN oin 

112. 


150 

RETURN 

MIN 0112' 

m. 



E'ln 

MIN 0113 
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SUBBOUTINE M*TV|C(*»i»C) 


1. 

2 . 

3 . C 

«. C THT 8 80 UTINC 8 fT 8 C«A TIHC 8 B 

5 - C C I R MU 8 T NOT Bt lOCNTICAl. VARIABLES 

t! 00 I HI , 3 

«• t cm«An,i)*ii(iHAn,?»*»( 2 )Mn, 5 )«i( 3 ) 

9, RETURN 

10, END 


1. SUBROUTINE M*TS(A,H,n 

2. RE»L A(S»T),H(3»3),C(3,8) 

3. C 

«. C THIS BOUTINE SET C«* TIMES H 

C r MUST NOT BE inrNTirAL ’«ITM either factor 

A, C 

7. on 1 in, I 

«r DO I .1*1,3 

< cf f,n«n.),n*ft(i,i)AA(j,?)*B(?,i) 

10. M ♦AIJ,T)*H(l,n 

It, RETURN 

12 . END 


SUBROUTINE MATMAT(A,B,C,L,M,N) 

rjN) - 


3 , C 

•8-;- - C THIS ROUTINE SETS C«A TIMES B “ * *' ' *• *‘“ 

5 . C C MUST NOT BE IO'ENTTCAL UlTH EITHER FACTOR 

h, -IT -- -- - - - ■ ■ 

7 , 00 I I*t,L 


— (j- ^ 

9. SAVE-0, 0 

to. DO 2 K* 1 ,M 

11. 2 SIVE»SAVE>A(1,K)*B(K,J) 

12 . -1 t(T,J)* 8 AVE 

IJ. RETURN 

i<r; * ' • END 
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# 


t 


1. SUPHOyTJNI 

2. J«KA,«K"#U) 

3. niMENSION 

U, HlLflC6,fc,6),AKAO(J#2),AKSnn,2),AlfA(5,2,6),A>'IUJ,2r6), 

5. 2BKACS»?»6)»l!Kpn,2,6),U(5,6) 

6. 00 JOO 

7. 00 SOO U>l#2 

8. 

•. AKAn(I,UJ«0. 

10. AK80(I»(..I)>0. 

11. 00 100 L«l»0 

12. AAAOn»W)«AKAP(M.J)*U(I.L)*AI AO(l.,J) 

13. lOO AK80(!,l.J)«AK80(I,l.J)*U{|,U*Al.H0(U»J) 

M. 00 ?0n K«l,fc 

II. AKA(I,LJ,K)1II, 

!£,. flHA(I.U»t^)«0. 

IT. AKH(I,lJ,K)«0. 

18. 8KBa,lJ.K)*ll. 

10. 00 200 l»l»6 

20. AKA(|,U,K)«*kA(THJ»tO*U{I,L)*A|.Aa»J*K) 

21. BKA(!,LJ»K)i«KAn.l.J,KWlin,L)*BLA(L,J»K) 

22. AMBd,! J,K)«AABn,U»KUU(I,U*A|.8a,J,K) 

2J. 200 BKR(I,|,J,K)iBKf.tr,U#K)*ll(T,U*BtHa,,a,K) 

20. JOO CONTTNIir 

25. f^ETUR'J 

26. END 




I. SUMOUTTNF 7«ERO(AKAO#AKAO«*Kt,Ai<B*RK«»nKf|.«(.PNA,PHl>CI.«* 

2* IPAMA»PAM2,AM) 

J. c 

A. C AM ■ AfROOVNAMIC TBPOUf 

5. C PIMA ■ PAPTIAL or AM WITH PiSPECT TP ALPHA 

fc. C PIMP ■ PAPTIAL or AM PITH PC8»ECT TO PHI 

7. C 

P. OIMCNSIOH AMf5),PAMA(S),PAMP(S),rxX(S,2),GXX(S»2),rKrP(S,2), 

9 , tCJOIP(5,?T,IKin(5.2)»AKHO(5.2)»IPI(5»2,P>.AKB(J*2»P)»IIKAn,2,B), 

10. 2BKR(S,2«6) 

11. PO UOO I«l»5 

12. AM(i)«f,, 

n. PAMAcn«o. 

la, pamp(ti*o. 

15. on 500 J«l»2 

IP. Finm,J)«AKA0(T.»I)/2, 

17. GIWU, !)«IKHn(I,.n/?, 

IP, FX*Pf!,J)»0. 

19. C*MP(I,J)«0. 

20. l.J«5*J/2 

21. XF«FL0AT(LJ) 

2?, rAL»CnpfXF*Al PHAI 

23. SAL«SIN(X7*ILPMA) 

29. pn 200 

25. TF«FLOAT(tO 

2P, CPH*C05(VF*PHI) 

27. SPH»3IN(YF*PHn 

2P. FXXPn,JI«FXXPf|,J)-Vr«UKA(I,J,M|*8PH-RKA((,J,»0*CPM) 

29. GXXP<1,J)»GXXP{I, n-Vr*( AKB( I,, >,X)«»8PH-PKB( I. J.KUCPH) 

50. FXX(l,.l)«rXXfI, J)#IKAn,J,K)«CPHlBKA{I,J,K)*SPH 

31 . 200 r.XXa..l)«t5XW(I..f)fAKP(t. J,K)*CPH«f}KPn,J,K)<i8PM 

52. PAMA(n*PIMA(n-XF*(rxY(T,J)*8IL-r.XX(I»J)*CAU 

33. PIMP(T).PAMpr ntFXXP(T.J)»CAURXXPn,J)*SAL 

30. 300 Awn),AMn)tF)(M(!,J)*CILir.XXfT,J>*8AL 

15. 900 CONTINUF 

50. 00 500 I«l,5 

37. PAMP<l)«aA<iPA»'Pn) 

5P. PAMA(I)*CLI*PAMA(n 

39. 500 AM(l)*Cl,A»AMn) 

90. BETIIHN 

91. e»JD 
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1 . 

2 . 

5 . 

tt. 

5. 

6 . 
T. 
«. 

10 . 

n. 

13 . 

i«. 

15. 

16. 
IT. 
Ift. 

19 . 
^0. 
21. 
2 ?. 
23. 

20. 
2‘.. 
26. 


sim»oijTiMf ns(*» wi. t, rvi, !vi» w» nm, ikm) 

C ^ISP*CK PROCrOUor FC’R all FICCNVALUES and eiGENVrCTOWS OF A RFAL 
C RENEWAL matr’X 
C A hEAL tNPilT matrix 

C NB REAL PART OF EIGENVALUES STORED IN A REAL vrCTOR 

C MI JMARINAPV PART OF EIGENVALUES STORED IN A REAL VECTOR 

C I MATRIX OF COMPLEX EIGENVECTORS STORED IN EISPACK 

C PACFEO FORM |N A REAL MATRIX 

C FVl WOP.KINC. VECTOR (REAL) 

C IVI rtORKINC VECTOR (INTEGER) 

C N NUMBER OF ROV*S AND COLS IN A 

C NM DIMENSIONFO col LENGTH OF A AND Z 

f lEPP nsPACK ERROR FLAG .NE. 0 AN ERROR MAS "EfN DFIECTED 

DlMENilOU A(UM, N), ?(NM, w), *«Rn), Pl(l), FVKI), IV|(1) 

CALL BfllANCrNM, N, Af ISI» IS2, FVl) 

CALL FIMHFS(NM, N, TSl , IS?» ». IVI) 

CALL ELTRAN(NM, N, TSI» IS2» A, (VI, Z) 

CALL HOR?(NM, N, ISl, IS?, A, HR, WI, Z, lERR) 

IFHERR .NE. 0) GO TO 999R 

CALL RALRAK(MM, N, ISl, IS2, FVl, N, Z) 

RETURN 

0Q99 PRINT 12, lERR 

\? FORMATOHO, <(5H**«**EPR0R in EISPACK EIGENVALUE ROUTINE HOR2 , 

* 2X, 6HIERR «, 16 ) 

RETURN 
' ENf) 
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SUBROimNE fUL«NC(NM,N»Af LOW, IGH, SCALE) 

69215009 

2. 

C 


69215005 

). 


INTEGER I,J,K,LtM,N,JJ.NM,I6H,L0W,lEXC 

69215006 

<1. 


REAL A(Nm,N),SCALE<N) 

69215007 

S. 


REAL C,E>C,R,S,R2«RA0I1i 

6921500B 

6. 

C 

real ars 

69215009 

7, 


LOGICAL NOCONV 

69215010 

s. 

C 


69215011 

9, 

C 

THIS SURROUTINE is a translation of the ALGOL RROCEDURC BALANCE, 

69215012 

10. 

C 

NUH, M*TH. 293-30«(1969) ST PARLFTT AND RETNSCM, 

69215013 

U. 

C 

HANDBOOK FOR AUTO, COMP., VOL. T I-LINfAR ALGEBRA, 315-326(1971). 

69215014 

12. 
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69215015 

13. 

C 

THIS SI'BROUTINF balances A REAL MATRIX AND ISOLATES 

69215016 

H. 

c 

EIGENVALUfS WHENEVER POSSIBLE, 

69215017 

15. 
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6921501B 

16. 

c 

ON INPUT- 

69215019 

17. 

c 


69215020 

18. 

c 

NM must be set to the row dimension of two-dimensional 

69215021 

19. 

c 

array PAPAMFTERS as oeclarfd in the calling program 

6921502< 

20. 

c 

DIMENSION STATEMENT, 

69215023 

21. 

c 


(■9215029 

22. 

c 

N IS THF OROFP OF THE MATRIX, 

69215025 

25. 

c 


69215026 

29. 

c 

A contains THE INPUT MATRIX TO BE BALANCED. 

69215027 

25. 

c 


6921502B 

26. 

c 

DN OUTPUT- , 

69215029 

27. 

c 


69215030 

2". 

c 

A contains THF BALANCED »*ATRIX, 

69215031 

29. 

c 


69215032 

30. 

c 

LO/» AND ir.H ape TWO INTEGERS SUCH THAT A( 1 ,J) 

69215033 

31. 

c 

!8 EQUAL TO ZERO IF 

69215034 

32. 

c 

(1) I IS GREATER THAN J AND 

69215035 

33. 

c 

(2) J«1 LOW-1 09 1«IGH*1,.,.,N, 

69215036 

39. 

c 


69215037 

35, 

c 

SCALE CONTAINS INFORMATION OFTFPVTNING THf 

69215038 

36. 

c 

permutations and SOALTNG FACTORS IISFO. 

69215039 

37. 

c 


69215040 

38. 

c 

SUPPOSE that the PRTNCIPAL submatrix IN ROWS LOW THROUGH IGH 

69215041 

39. 

c 

HAS BEEN PAI.ANCFO, THAT P(J) DENOTES THE INDEX INTERCHANGED 

69215042 

90. 

c 

WITH .1 DUPING THE PtPMUTATTON STEP, AND THAT THE ELEMENTS 

69215043 

91. 

c 

OF THE DIAGONAL MATRIX USED ARE DEMOTED BY 0(1, J). THFN 

69215094 

92. 

c 

scAiF(j) « P(.n, rop J * i,...,LOw-i 

69215045 

93. 

c 

S J « LOW,...,IGH 

69215046 

99. 

c 

= P(J) J « ighti,.,.,n. 

69215047 

95. 

c 

THE ORIfR IN HHICH THE TNTFRCHANGES APE MADE IS N TO IGHTl, 

6921504B 

96. 

c 

THEN 1 TO LClH-t. 

69215049 

97. 

c 


69215050 

98. 

c 

,'!OTE THAT 1 IS RETURNED FOR IGH IF IGH IS ZERO FORMALLY, 

69215051 

99. 

c 

. 

69215052 

50. 

c 

the ALGOL PROCEDURE FXC CONTAINED IN BALANCE ARREARS IN 

69215053 

51. 

r. 

RALANC in line, (NDTF that the ALGOL ROLES OF lOpNTlFlERS 

69215054 

52. 

c 

K,L HAVE BEEN REVERSED.) 

69215055 


»s. 

C 







69215056 

»•. 

c 


oufSTinwj am cw^'ents should 

BE directed to ». S. 

CARBOW, 

69215057 

55. 

c 


*PELieO M*THfM*TICS DIVISION, 

ARCONNE NATIONAL LABORATORY 

6921505B 

56. 

c 







69215059 

57. 

c 



dependent 

RARAHETER 

SPfCTfYING 

69215062 

56. 
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THE BASE or THE 

HACMINC FLOATINC ROINT 

RERRFSCNTATIDN. 

69215065 

5*. 

c 







69215069 

60. 

c 







69215065 

61. 



RADIX a 2. 





69215066 

62. 

c 
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65. 



82 a radix * RADIX 





69215068 

60. 



K a 1 





69215069 

65. 



L a N 





69215070 

66. 



on Tfi 100 





69215071 

67. 

c 




69215072 

66. 

c 


COLUMN EXCHANGE 




69215073 

66, 


20 

SCALEtM) « .1 





69215079 

70. 



IK (J .ED. M) GO TO 50 





69215075 

7!. 

c 







69215076 

72. 



DO 30 I a 1 , L 





69215077 

75. 



r a A(I,J) 





69215078 

7<<. 



A(1,J) * A(I,M) 





69215079 

75. 



A(I#M) * F 





6R215060 

76. 


30 

CONTIMlF 





69215081 

77. 

c 







69215062 

76. 



no 60 1 a K, N 





69215083 

76. 



F « A(j,n 





69215089 

SO. 



A(j,n « A(M,n 





69215085 

81. 



A(M,n a F 





69215086 

82. 


60 

CONTINUE 





6R215087 

83. 

r 
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66, 


50 

GO TO (80,130), TCXC 





69215089 

85. 

c 


*••**«**•* SEARCH FOR RONS 

isolating an 

EIGENVALUE 

69215090 

86. 

c 


ANP PUSH THEN DOWN 




69215091 

87. 


80 

IF (L .EO. 1) GO TO 280 





69215092 

88. 



L a L - 1 





69215093 

86. 
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********** FOR JaL STEP -1 

UNTIL 1 00 -« 


* 

69215099 

60. 


100 

DO 120 JJ a 1, L 





69215095 

61. 



J a L ♦ 1 - JJ 





69215096 

62. 

c 







69215097 

65. 



00 1 10 I s j , L 





69215098 

66, 



IF (I .FO. J) GO TP 

no 




69215099 

65. 



IF (A(J,I) ,NE. n.O) 

GO 

TO 120 



69215100 

66. 


110 

CONTINUE 





69215101 

67. 

c 







69215102 

68. 



M a L 





69215103 

66. 



IF.XC a I 





69215109 

100, 



GO TO 20 





69215105 

101. 


1?0 

CONTINUE 





69215106 

102. 

c 







69215107 

105. 



GO TO 160 





69215108 

106. 

c 



ISOLATING 

AN EIGENVALUE 

69215109 
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105, 

C 


UNO PUSH TMFH L€F7 

• •• 

• ••« 


1 

6921SI10 

106. 


130 

K ■ K 6 1 





69215111 

107. 

c 







692IS1I2 

I0«. 


too 

no 170 J ■ K, L 





69215113 

109. 

c 







69215116 

no. 



DO ISO I ■ K, L 





69215115 
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TF (I .EU, J) GO TO 150 





69215116 

112. 



IF (A(I,J) .fiE, 0,0) GO 

TO 

170 



69215117 

115. 


150 

CONTINUE 





69215116 

n<i. 

c 







69215119 

ns. 



M ■ K 





69215120 

116. 



lEXC » 2 





69215121 

117. 



GO TO ?« 





69215122 

ns. 


170 

CONTINIIF 





69215123 

119. 

c 
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**•* 69215126 

120. 



00 ISO I « K, L 





69215,25 

121. 


ISO 

SCALEH) « 1.0 





69215126 

122. 

c 


««#•***•*« TTF9ATIVE ino® FOR 


69215127 

125. 


190 

NOCONV « .PAUSE. 





69215128 

120. 

c 







69215129 

125. 



no 270 1 ■ 6, L 





69215130 

126. 



c « 0,0 





69215131 

127. 



R ■ (>.0 





69215132 

12«. 

c 







69215133 

129. 



DO 200 J m K, 1 





69215136 

150. 



TF (J-.FO, T) GO TO 200 





69215135 

151 , 



C * C ♦ A8S(A(J,I)) 





69215156 

152. 



P ■ P ♦ AOS(A(I,.in 





69215157 

155. 


206 

CONTINIIF 





69215158 

150. 

r. 







69215159 

135. 



C, • N / RAnjX 





69215190 

136. 



F ■ 1,0 





69215161 

137. 



S ■ C ♦ R 





69215162 

15S. 


210 

IF (C ,r,E, r.) GO TO 220 





69215165 

139. 



Fir* RADIX 





69215196 

100. 



C « C • B2 





69215165 

1«1. 



GO TO 210 





69215166 

102. 


220 

G ■ R * RADIX 





69215167 

105. 


250 

IF (C .UT. G) GO TO ?00 





69215168 

100. 



F * F / RAOTX 





69215169 

105. 



C « C / R2 





69215150 

lOfc. 



GO TO 250 





69215151 

107. 

c 



*•* 




69215152 

ion. 


20 0 

IF ((C ♦ 9W f .GE. 0.95 * 

S) 

GO 

TO 

270 

69215155 

109. 



G s 1.0 / F 





69215156 

150. 



SCALE(I) ■ SCALF(I) • F 





692I5I5S 

151. 



NOCONV a .TRUE. 





69215156 

152. 

c 







69215157 

153. 



DO 250 J ■ K, N 





69215158 

ISO. 


250 

A(T,J) « A(I,J) * G 





69215159 

155. 

c 







69215160 

156. 



on 260 J * n u 





692I5I6I 

157. 


260 

A(J,I) e A(J,!1 * r 





69215162 

158. 

c 







69215165 

159. 


270 

CONTINUf 





69215166 

160. 

c 







69215165 

161 . 



IF (NOCONV) GO TO 190 





69215166 

162. 

C 







69215167 

165, 


260 

LOW « K 





69215168 

160. 



IGH s U 





69215169 

165. 



RETURN 





69215170 

1 66. 

c 


.********* last CARO np RALANC 

*«#*#*«•** 

69215171 

167. 



Fv.n 





69215172 
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SUBROUTIHR ELMHfS(NM,N,LOW,ICH,*,I»<T) 

INTEOER T,J,m,N,LA»NM,1RH,KP1,LOW*MMJ ,HP| 

REAL AfNM,M) 

REAL X»y 
REAL ARS 
INTEGER INTIIRM) 

THIS SUBROUTINE IS A TRANSLATION OF THE ALGOL PROCFOHRE ELMMES# 
MUM, m«th„ 1?, 50R-368(1R<.«) BY MIRTIN AND WILKINSON, 

HANDBOOK FOR AUTO. COMP,, VOL . I I-L INE AR ALGEBRA# SSR-35B(IR7I). 

GIVEN A PEAL OENFRAL MatRIY, THIS SUBROUTINE 
REDUCES A SURMATRIX SITUATED IN ROWS AND COLUMNS 
LOW THROUGH IGH TO UPPER HESSENBFRG FORM «Y 
stabilized FLf'MFNTARY SIMILARITY TRANSFORMATIONS. 

ON INPUT- 

NM MUST HE SET TO THF POrt DIMENSION OF TWO-OIMENSTONAL 
AMRAr PAPAMFTFRS AS DECLARfO IN THF CALLING PROGRAM 
OIMfNSlON STATFMENT# 

N IS THE OROFR OF THF MATRIX, 

LOW ANr tRH ARF INTFGF^S DFTER^'INFO BY THf BALANCING 
S'iDPOUTI'iF PALANC. IP RALANC HAS NOT RFFM USED, 

SET town, TGHrN, 

A contains the INPUT matrix. 

ON OUTPliT- 

A contains THF HESSFHHFRG MATRIX. THE MULTIPLIERS 
WHICH HERE USFD TN THE RFOUCTIOM ARE STORED TN THE 
REMAINING TRIANGLE UNDER THE HESSENBFRG MATRIX, 

INT CONTAINS information 0^ THE RO>'S AND COLUMNS 
interchanged in THF REDUCTION, 
only ELEMENTS LOW THROUGH IGH ARE USED. 

OUESTIONS and comments should PE DIRECTED TO H. 8. GARrOW, 
APPLIED mathematics DlVISjnf, APGONNE NATIONAL LABORATORY 


LA s IGH - 1 
KP 1 s LOW ♦ 1 

IF (LA .LT. KPt ) GO TO ? 0 n 

DO ISO M ■ KP1 , la 
MMJ ■ H . I 
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ss. 



t ■ O.A 

7S21S056 

5<». 



I ■ M 

75215057 

55. 

c 



71215658 

56. 



DO IBO J ■ M, ICH 

75215059 

57. 



ir (6BS(*(J»MHn) ,LF. *«S{»<)) GO TO mo 

75215060 

SB. 



t m 

75215061 

59, 



I ■ .1 

75215062 

60. 


100 

CONTINUE 

73215065 

6i. 

c 



73215069 

62. 



INT(“) ■ I 

75215065 

65. 



IF (I .EO, M) GO TO no 

75215066 

69. 

c 

i 


75215067 

65. 



on 1 10 J * mmi . M 

75215068 

<>6. 



V * 

75215069 

67. 



A(T,J) * A(W,J) 

73215070 

6fl. 



« V 

73215071 

69. 


vio 

CONTTMIt 

75215072 

70. 

c 



75215075 

71 . 



nn ipo j X 1 , ir.H 

75215079 

7?. 



V « A(,I,T) 

75215075 

7*i. 



A(J,n « A(J,M) 

75215076 

7/J. 



« V 

73215077 

75. 


l?0 

CnNTINUF 

7521507" 

76. 

c 

i 


75215079 

7?. 


no 

IF (X .Fo. 0,0) GO rn mo 

73215080 

7«. 



MPJ 1 M ♦ 1 

75215081 

7«. 

r 



75215082 

HO, 



r>0 I6U J • '*Pl, ir.H 

75215085 

«l . 



V = An.MMi ) 

73215089 

H2. 



IF (Y .Eo. 0 . 0 ) r.n Tn mo 

75215085 

«3. 



Y S Y / X 

73215086 




1 ) e V 

75215087 

65. 

c 



75215088 

B6. 



no 100 J r M, 

75215089 

A7. 


too 

A(I,J) X AH.J) - Y * A(M,J) 

75215090 


c 



73215091 

HQ. 



00 150 J S 1» IGH 

73215092 

90. 


150 

A(.T,M) ■ A(J,M) ♦ Y * A(J,n 

75215095 

91. 

c 



73215099 

9?. 


160 

CONTINUE 

15215095 

93. 

c 



73215096 

90. 


1 HO 

CONTTNIIF 

73215097 

95. 

c 



73215098 

96. 


?no 

RETURN 

75215099 

97. 

r. 

i 


73215100 

9H. 



ENO 

75215101 


s 

i 



t. 

2 . 

S. 

«. 

S. 

6 * 

7. 

10 . 

n. 

1 ?. 

13. 

10 . 

IS. 

u. 

17. 

»H. 

l*». 

20 . 

21 . 

2 ?. 

23. 

20 . 

2*i. 

26. 

27. 

20 . 

20 . 

30. 

31. 

32. 

33. 
30. 

35. 

36. 

37. 
30 . 
39. 
00 . 
01 . 
02 . 
03. 
00 . 

05. 

06. 

07. 

08. 
09. 

50. 

51. 

52 . 


C 


C 

C 

C 

C 

C 

r 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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c 

r 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

r. 

c 

c 

c 

c 

c 
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SUWOUTINC RtTH0NCNM.N,L0W,|CH,4,INT,Z) 

INTC6ER J^N,KL,MM«MP,Nr4.IGN.L01«,MPl 
Rr*l A(NM,IGH),Z(NM,N) 

INTCCCR TNT(If.M) 

THIS SI'BROUTiNr !S * TRAMSLATinN OF THE OLGOl FHOCEDUFC €LMTf»«N8* 
NUM, MATH. |6» 181-200(1970) BY PETERS AND WILKINSON. 

HANDBOOK FOR AUTO. COMP.. VOL . I I-L INEAR ALGEBRA, S7?-SfS(|97n. 

THIS SIIRROUTINF ACCUMULATFS THE STARILIZFO ELEMENTART 

sihraritv transformations used in the reduction of a 

PEAL GENERAL MATRIX TO UPPER HESSFNHERG FORM BY ELMHFS. 

ON TNPUT- 

NM WUST BE SFT TO THE ROh ni«FN8I0N OF TnO-OIMENSIONAL 

array para»'ftehs as ofclarfo in the calling program 

DIMENSION STATCMFNT, 

N IS THE ORf^EP up THE MATRIX, 

LON ANO IGH ARf integers nETERMTNEO BY THE BALANCING 
SUBROUTINE BALANC. IF BALANC HAS NOT BEEN USED, 

SET L0N=1, IGHsN, 

A‘CONTAI^S THE MULTIPLIERS WHICH *«FRE USFO IN THE 
REDUCTION BY ElMHfS IN ITS LOWER TRIANGLE 
BELOW THE SUBOIAGONAl, 

INT CONTAINS INFPRMATlfN ON THF ROWS AND COLUMNS 
IF rFRCFiAr GlI) IN THE FFOUCTION BY ELMHES. 

ONLY ELEMENTS LOW THROUGH IGH ARF USED. 

ON OUTPUT- 

Z CONTAINS THF transformation MATRIX PPOOUCFP IN THE 
reduction by FLMHFS. 

OUESTIONS AND COmmfNTS SHOULD BF DIPECTFD TO R. S. GAPROW, 

APPLIED MATHEMATICS DIVISION, ARGONNF NATIONAL LABORATORY 


***«•*•*•* initialize 7 TO IDENTITY MATRIX 
DO »0 I » 1 , N 

DO e>n j B i, M 

60 . zn,j) = 0.0 

zd.ii = 1.0 
80 continue 
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51. 

C 



20215056 




KL ■ If.H - LOW - i 

20215057 

5S, 



IF (KL .LT. t) fiO TO ?00 

20215050 

56. 

C 



aaa* 20215050 

5T. 



DO HO a 1, KL 

20215060 

5«. 



Mp a IGM • MH 

20215061 

5«». 



MPt a Mp ♦ 1 

20215062 

60. 

c 



20215061 

6|, 



DO 1(10 I a MPI, !GH 

20215060 

62. 


1 10 

2(1, MP) a *(I,MP-1) 

20215065 

63. 

c 



20215066 

6«. 



t « IfJT(MP) 

202)5067 

t>5. 



IF (I .FO. MP) GO TO Mifl 

20215060 

66. 

c 



20215069 

67. 



DO ISO J a MP, IGH 

20215070 

6". 



2(MP,J) a 2(1, J) 

20215071 

6‘». 



7(1, J) ■ 0.0 

20215072 

70. 


no 

CONTIMIF 

20215073 

71 . 

c 



20215070 

72. 



Z(1,“P) • 1.0 

20215075 

73. 



CONTI Ml IF 

20215076 

70. 

c 



20215077 

1 %. 


?Ofi 

RFTlJWK 

20215070 

76. 

r 



20215079 

77. 



END 

20215000 


SUBROIITINF HQR2(N»«fN,lOW.ICH>H.WR,MT,{»ICBR) 

TNTEGER I, J,K,L,M,n,EN,1I,JJ,LL,MM,N*,NH,NN, 

X IGH#IT8,LOW,»'P2,ENM2»fEBH 

RE»L H(NM,M),WH(N)»WI(N),Z(NM,N) 

«E»L P,Q,»,S,T,W,X,V,B*,S*,VI,VR,Z2,N0RM,M*CHEP 

REAL SQHT.ABSfSIGH 

IfJTEGER M|N0 

LOGICAL NOTI AS 

CnPPLtX Z5 

COMPLEX CMPLX 

heal T3(?) 

EQUIVALENCE (Z3iT3(U) 

THIS SUHHOtlTINF IS A TRANSLATION OP THE ALGOL PPOCEPURp HOR2» 
N')M, MATH. 16, l8l-20A(l97(n BY PETERS AND WILKINSON, 

HAHOBOOK for auto, comp,, VOL.II-LINEAR algebra, 372-595(1971). 

THIS StlPROliIINf FINDS THF EIGENVALUES AMD EIGENVECTORS 
OF A real liRPEP HESSENBERG MATRIX BY THE QR METHOD, THE 
FIGENVECTOHS OF A REAL GENERAL MATRIX CAN ALSO BE FOUND 
IF ELMHFS AND ELTRAN OP ORTHES AND ORTRAH HAVE 
BEEN USEO TO REDUCE THIS CFNERAL MATRIX TO HFSSE.NHEPG FORM 
AMO TO ACCUMULATE THE SIMILARITY TRANSFORMATIONS, 

ON INPUT- 

MM MUST PF SFT TO THE RPt" DIMENSION OF TWO-DIMENSIONAL 
array PARAMFTERS as DFn.ARFD IN THE CALLING PROGRAM 
OIMFNSION statement, 

N IS THF ORDER OF THE MATRIX, 

LOrt and IGH are irnFGERS DETERMINED BY THE BALANCTNC, 
SIJOROUTINE PALANC,. IF PALANC HAS NOT RFfN USED, 

SET inw«l, IGHsM, 

H CONTAINS THF UPPER HFSSENBERG MATRIX, 

Z contains THE THAmsfORMATION matrix PRODUCED BY ELTRAN 
AFTFP THE REDUCTION BY FLMHFS, OR BY ORTRAN AFTER THf 
RFfHlfTION BY npTMES, IF PERFORmfD, IF THE EIGENVECTORS 
OF THE H( SSENPERG MATRIX ARE OpSIREO, Z MUST CONTATN THE 
IDENTITY MATRIX, 

PM OIITPUT- 

H HAS BFEN DESTROYED, 

WR AND W1 contain THE REAL ANO IMAGINARY PARTS, 

respectively, DF THF EIGENVALUES. THE EIGENVALUES 
APE IINOROFRED FXCEPT THAT COMPLEX CONJUGATE PAIRS 
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C 


OF VA1UF.8 6PPF4H CONSICUTIVFLV WITH THF ElCCNVALUt 

BT2I5956 


c 


H6VIN6 THF positive IMACIN6PV P6RT FliST. IF AN 

B72I505T 

55. 

c 


EOPOP FXn IS MAOE» THE ETCfNVALUCS SHOULD HE COMCCT 

87215058 

56. 

r. 


FOR INDlCfS IERPM«...«Nr 

87215059 

57. 

c 



87215060 

5«. 

c 


Z contains THF REAL AND IHAOINART PARTS OF THE EICCNVECTORS. 

87215061 

59. 

c 


IF THE l-TM ElfifNVALUE IS REAL* THE l»TH COLUMN OF Z 

67215062 

60. 

c 


contains its eigenvector, if the i-th eigenvalue is complex 

87215063 

6t. 

c 


WTTH POSITIVE IMACINAPY PART, THIJ |-TH AND (I6I)«TH 

67215069 

62. 

c 


COLUMNS OF Z CONTAIN THE REAL AND IMAGINARY PARTS OF ITS 

87215065 

63. 

c 


CIKENVECTOP. THE EIGENVECTORS ARE UNNORMALIZEO, IF AN 

87215066 

6«. 

c 


ERROR EXIT IS MAOF, NONE OF THE EIGENVECTORS MAS REEN FOUND, 

87215067 

65. 

c 



87215068 

66. 

c 


lERR IS SET TO 

87215069 

67. 

c 


ZERO FOR NORMAL RETURN, 

87215070 

66. 

c 


J IF THE J-TH EIGENVALUE HAS NOT BEEN 

87215071 

69. 

c 


DETCRMINEO AFTER 30 ITERATIONS. 

87215072 

70. 

c 



87215073 

7t, 

c 


arithmetic IS REAL EXCEPT FOP THF REPLACEMENT OF THE ALGOL 

87215079 

7?. 

c 


PROCEDURE CniV BY COMPIEX DIVISION, 

87215075 

73. 

c 



87215076 

7«. 

c 


DUFSTIDNS AND COMMENTS SHOULO 8F OIRECTEO TO H, 8, GARpnw, 

87215077 

75. 

c 


APPLIED mathematics DIVISIOM, ARGONNF NATIONAL LABORATORY 

87215078 

76. 

c 



87215079 

77 
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^ . .. . ^ ........ .... .... r ■ 1 ■ ^ ■ I" r ■ 


78. 

c 



87215081 

79. 

c 



87215082 

«o. 

c 


THE RELATIVE PRECISION Of FLOATING POINT ARITHMETIC. 

87215083 

*1. 

c 



87215089 

82. 

c 



87215085 

83. 

c 


sigma ^ SIGMA 5 SIGMA 5 S!G*-.45 SIGMa 5 SIGMA 5 SIGMA N SIGMA 5 

• • « A 

fiu. 



macmEP * 2.0«*(-52) 


85. 

c 


SIGMA 5 Sir.MA 5 SIGMA 5 SIGMA5 SIGMA 5 SIGMA 5 SIGMA 5 SIGMA 5 

• A A A 

86. 

c 



87215087 

"7. 



IfRH * 0 

87215088 

88. 

c 



87215089 

89. 



DO 50 I * 1, N 

87215C90 

90. 



IE (I .GF, LOW .ANO. 1 .LE. IGH) GO TO 50 

87215091 

91. 



WR(I) * H(I,1) 

87215092 

9?. 



wun « o.n 

87215093 

93. 


50 

COMTINIlt 

87215099 

OU. 

c 



87215095 

95. 



EN = IGH 

87215096 

96. 



T » 0,0 

87215097 

97. 

c 



87215098 

98. 


60 

IF (EN .LT. LOV.) GO TO 390 

87215099 

90. 



ITS * 0 

87215100 

100. 



NA * EM - 1 

8T21510I 

101 . 



ENM3 » NA - 1 

87215102 

102. 

c 



87215103 

103. 

c 



67215109 

inu. 


70 

DO 80 LL « LOW, EN 

87215105 
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las. 



L ■ EN ♦ LOW • LL 

872ISI06 

IH. 



IF (L ,EO, LOW) 60 TO 100 

BT2I1107 

l(?i, 



ir t688{H(L»L-m .LE, MICHfP * («oa(H(L»l ,L*t ) ) 

872lffOB 



X ♦ ABSIHIUL)))) 00 TO }O0 

B721510B 

»')9, 


60 

CONTINUE 

B7215I10 

tU^i 

C 



B72I5UI 

m= 


too 

X ■ H(fN,EN) 

87P15H2 

H2. 



IF (U .ro. EN) CO TO 270 

B7215113 

113. 



Y M 

67215119 

119. 



M ■ H(ENfNA) * H(N4,FN) 

87215115 

113, 



IF <L .10, N4) r.O TO 280 

B7215I16 

116. 



IF (ITS ,EO. 50) GO TO |000 

87215117 

117. 



IF (ITS ,NF, 10 .4N0, ITS .6E, 20) 60 TO 150 

87215118 

tIA. 

C 



67215119 

ll'?. 



T » T Y X 

§7215120 

120. 

C 



87215121 

121. 



00 120 I « LOW, EN 

67215122 

122. 


120 

H(I,I) ■ H(f,|) - X 

67215125 

J25. 

C 



87215129 

129. 



S ■ ABS(H(£N,N«n ♦ JlHS(«(f'A,FNM2)) 

87215125 

I2S. 



X • 0,75 * S 

87215126 

126. 



Y * X 

67215127 

127. 



6 * -0.9375 * S * S 

87215126 

126, 


130 

ITS * ITS ♦ 1 

67215129 

129. 

r 



87215130 

130. 

c 


suB-niAGONiL elements. 

67215131 

151, 

c 



87215132 

132. 



on 190 MM « L, ENM2 

87215133 

133. 



M « ENM2 ♦ L - MM 

87215139 

139. 



Z7 * H(M,M) 

67215135 

153. 



9 « X .. 27 

67215136 

136. 



s « y - 7Z 

6.7215137 

137. 



p S (R * s - W) / ♦ H(M,Mt!) 

67215136 

136. 



0 « H(M*i ) . 7Z . f . g 

67215139 

139. 



p « 

87215190 

190. 



S * AHS(P) ♦ AHS(O) ♦ SRS(R) 

67215191 

191. 



P » P / s 

67215192 

192. 



0 s fJ / S 

67215193 

193. 



t? » M / S 

87215199 

199. 



IE f-/' .EO. U GO TP 150 

67215195 

195. 



IF * (ABSfOl ♦ AHS(P)) ,LE. ^acHEP * AflSfP) 

67215196 

196. 



X * (Af»S(H(M-t ,M-in ♦ APSfZZ) ♦ ABS(H(»'+l,M*m)) GO TO 150 

67215197 

197. 


1/10 

COMTIVUt 

07215196 

196. 

c 



67215199 

199. 


150 

MP2 = M ♦ J 

67215150 

150. 

c 



6721515! 

151 . 



pn 16(1 1 a Mp?, fK 

67215152 

152. 



H( 1,1-2) « 0.0 

87215155 

153. 



• IE (I .EG. (-tPai GO TO 160 

67215159 

159. 



H(l,I-3) « 0.0 

67215155 

155. 


160 

CONTINUE 

67215156 

156. 

c 


•**««•**** POUBLE OR STEP INVOLVING ROWS L TP FN AND 

67215157 
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157. 

C 


190, 


DO 260 K ■ N« 

154, 


NOTIAS ■ K .NC. NA 

160, 


IF (K .rO. CO TO 170 

161. 


p I H(K,0-U 

162. 


0 ■ H(KTUH-l) 

1*1. 


R ■ 0.0 

164. 


IF (NDTLAS) R ■ H(K62,K.t) 

165. 


X a 4HS(R) ♦ ARS(Q) ♦ A00(R) 

166. 


IF (X .FO, 0,0) GO TO 260 

167, 


P m P / n 

1*0, 


0 a 0 / K 

164. 


R a R / X 

170, 

170 

8 a 81RN(SORTfR*R*Oa(»*R»R),R) 

171. 


IF (0 .CO. **) 00 TO 100 

172. 


H(K,K-n a .$ * X 

175. 


GO TO 190 

174. 

100 

IF a ,NC, M) h(K,K-1) a -H(k,K-l) 

179. 

140 

P a P ♦ s 

176. 


X a P / 5 

177. 


V a W / 8 

170, 


U • P / $ 

174. 


0 a 9 / P 

100. 


R a R / P 

101. 

I"2, 

C 


DO 210 J a K, f 

105. 


P a H(0,a) ♦ f) • HfK^l.J) 

104. 


TP (.^'OT. NOTLASl GO TO 200 

109. 


P ■ P ♦ P * HTKT2, I) 

106. 


H(K*2,J) a H(Kt2,Jl - P • 72 

107, 

200 

H(K*l,,l) • H(Kf),J) - P * ¥ 

100. 


H(K,J) a - P • X 

169. 

210 

CONTINDF 

190. 

C 


141. 

142. 
145, 

r 

J a VINn(EN,X*5) 


DO 250 I • 1, J 

194, 


P ■ X * Hfl.K) 4 Y * M{I»K41) 

149. 


IF (.MPT, nnn.*S) GC TO 220 

146. 


P a P 4 77 4 H(1,K42) 

147. 


‘^(I,K421 a H(I,K42) - P • P 

140. 

220 

H(I,K4l) a H(1,K4|) - P * 0 

149. 


H(I.K) a H(I,K) - P 

200. 

250 

CPNTIMJF 

201. 

C 


202. 


DP 290 I a LDiT, IGH 

203. 


P a X 4 Zn,K) 4 Y * 2(1, K41) 

204. 


IF (.^OT, NflTLAS) GO TO 240 

200, 


P a P 4 27 4 7(1, H42) 

206. 


2(t,K42) a 7(1, K42) - p * R 

207. 

240 

7(1, K4I) a 2(1, K41) - P 4 0 

200. 


7(1, K) a 7(1, K) - P 


mifisi 

STtlSUf 

072191*0 

072191*1 

072191*2 

072191*1 

072191*4 

072191*9 

072191** 

072191*7 

072191*0 

072191*4 

07219170 

07219171 

07219172 

0T3I917J 

07219174 

07219179 

07219176 

07219177 

07219170 

07219174 

67219100 

07219161 

07219102 

072191OS 

07219104 

07219109 

07219106 

07219107 

07219100 

07219104 

67219140 

07215141 

07215142 

07215193 

07215144 

07215149 

07215146 

07215147 

07215190 

07215144 

67215200 

07215201 

07219202 

07215205 

07219204 

07215205 

07219206 

07219207 

07219200 

07215204 
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209. 


250 

CONTI DUf 

87215210 

210. 

C 



87215211 

2H. 


260 

CONTINUE 

87215212 

212. 

c 



87215213 

211. 



CO TO 70 

87215216 

210. 

c 



87215215 

215. 


270 

H(CN,ENJ • X ♦ T 

87215216 

2 1 6 • 



NIMtN) ■ H(fN.FN) 

87215217 

217. 



VnI(fN) ■ 0.0 

87215218 

210. 



EN « Hk 

87215219 

219. 



50 TO 60 

87215220 

220. 

c 



87215221 

221. 


?80 

0 ■ (V - X) / 2,0 

87215222 

222. 



0 * P • P ♦ X< 

87215223 

225 



u « sopTt6BS(nn 

87215226 

220. 



H(FN,|N) « X ♦ t 

87215225 

225. 



X » M(fN.EN) 

87215226 

220. 



H(NA,N6) « Y ♦ T 

87215227 

227. 



If U4 .IT. 0.0) on TO 320 

87215228 

22«. 

c 


*••***••*• PfXL 9A10 ***«*****• 

87215229 

2?o. 



ZZ * p f 9ir.N(ZZ#P) 

87215230 

250. 



wR(N») B X ♦ 77 

87215231 

231. 



WP(fN) « wM(NA) 

87215232 

233. 



IF (ZZ .NF. 0.0) 5MUF.N) » X - / ZZ 

872)5233 

233. 



NIfVA) s 0.C 

87215236 

230. 



’YI(FN) B 0,0 

87215235 

23S. 



X B H(FN.UA) 

87215236 

236. 



R ■ S0«T(X*X»7Z*77) 

87215237 

237. 



P B X / R 

87215238 

23". 



0 B Z7 / « 

87215259 

230. 

c 



872)5200 

200. 



00 200 J * NA, N 

87215261 

201. 



77 * H(MA,.n 

87215262 

202. 



HCNA.J) B 0 * ZZ f P * 

87215263 

203. 



HCEKmD « Q * - P * ZZ 

87215266 

200, 


200 

CONTINUE 

87215265 

20S. 

f 



87215266 

206. 



on 300 I s 1 , FA' 

87215267 

207, 



ZZ s Hd.AlA) 

87215268 

208. 



H(l,NA) * fj * 77 ♦ P * H(I,FN) 

8721526R 

200. 



Hn.EN) B D * H(I,FN) - 0 * 7Z 

87215250 

250. 


500 

CON’IMit 

87215251 

251. 

c 



87215252 

252. 



on 310 I e LOW, ICH 

87215253 

253, 



ZZ = Zn.NA) 

87215256 

254. 



Zn.AiA) B f) * 77 ♦ P • 7fI,EN) 

87215255 

255. 



Z(I,f'4) * 0 • 7(1, FN) - P * 77 

67215256 

256. 


310 

CONTINUE 

87215257 

257. 

c 



87215258 

250. 



no TO 330 

87215259 

250. 

c 



67215260 

260. 


320 

NR(NA) B y ♦ P 

87215261 
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261. 



M9fCN) ■ X ♦ 2 

072!52«2 

> 262. 



2!<N6) ■ 22 

07215265 

263. 



Wl(CN) ■ -22 

97215266 

26tt. 


550 

CN 

■ EN22 

97215265 

265. 



60 

TO 60 

97215266 

266, 

C 


tm 


97215267 

267. 

C 




97215269 

26«. 


3«0 

N06M a 0,0 

07215269 

? 269. 



K a 1 

97215270 

^ 270. 

C 




97215271 

271. 



no 

560 1 a 1, N 

97215272 

272. 

C 




97215275 

1 275. 




00 550 J a K, 

97215276 

f 27U, 


350 


MPRM a NOPM ♦ ABSfHd.J)) 

97215275 

' 275. 

C 




97215276 

276. 




K a I 

97215277 

277, 


5*jO 

CONTI'JUF 

97215279 

1 276. 

C 




97215279 

, 279. 



IF 

(NQHM .10. 0.0) GO Tn 1001 

97215290 

: ?«o. 

C 


• *< 


97215291 

251. 



no 

ftOO NN a 1, N 

B7215282 

252, 




a 1 ♦ 1 • TJfj 

87215293 

253. 




9 a nR(F'J) 

97215286 

2»U. 




0 a wl(H‘l 

87215295 

255. 




VA a F*j • 1 

87215296 

266. 




IF (n) 710. 600 , 800 

87215287 

257. 

C 


• *( 


87215288 

256. 


600 


M a ffi 

87215289 

259. 




H(tN,EN) « 1,0 

87215290 

?9n. 




IF (►iA .EG. 0) on TO 800 

87215291 

291. 

C 


* 44 


872l!i292 

292. 




nn 700 n * 1, NA 

87215293 

293. 




I * EN - IT 

87215296 

29«. 




n a H(I,n - P 

87215295 

295. 




9 a 

87215296 

296, 




IF (M ,r,T. MA) GO TO 620 

87215297 

297. 

c 




87215298 

298. 




00 610 J 5 V, NA 

87215299 

?99. 


610 


B « R ♦ H(I,J) * HIJ.EN) 

87215500 

300. 

c 




87215501 

301. 


620 


IF (WKI) .r,E. 0,0) r,n in 630 

87215302 

502. 




n a W 

87215305 

303. 




sap 

87215506 

30/J. 




GO TO 700 

87215505 

305. 


630 


M a I 

87215306 

30(S, 




IF (Wfim .^'F.. 0 . 0 ) nn TO 600 

87215307 

30 7. 




Taw 

87215508 

30R. 




IF fW .Ft). 0.0) T a maCHFP * MORM 


J09. 




H(I,F6) = -H / T 

'id 

310. 




GO TO 700 


Ml. 

c 



/J4i‘5512 

M2. 


600 


* = 11(1, 1 + U 

»7J;53jS 


38 


j 


MS. 



V * HUfl,!) 

87215314 

310. 



0 ■ (WR(n - p) 9 (6R(i) - P) ♦ wi(i) • wim 

87215315 

315. 



T ■ fX * S • ZZ • R) / 0 

87215316 

316. 



H(I,CN) 1 T 

87215317 

317. 



IF (APS(X) ,Lt, ABS(ZZ)) GO TO 650 

87215318 

318. 



Hdfl.EN) * (-R - W * T) / X 

87215319 

S19. 



GO TO 700 

87215320 

320. 


650 

H(I41,EN) ■ (-S - V t tw zz 

87215321 

321. 


700 

CONTINUE 

87215322 

322. 

C 



87215323 

323. 



00 TO "00 

87215324 

32«. 

c 



87215325 

325. 


710 

M 9 NA 

87215326 

326. 

c 



87215327 

327. 

c 



87215328 

328. 



IF {ABS(H(FN,NA)1 .LF. APS(H(NA,EN)>) go to 720 

87215329 

329, 



H(NA,NA) 10/ H(EN,MA) 

87215350 

330. 



PlNA.fN) ■ -(M(EN#tN) - P) / H(EN,MA) 

87215331 

331. 



GO TO ’30 

87215332 

332. 


7?0 

Z3 ■ C*'PLX(0.0,-H{NA,EM) / CMPLX(H(NA,NA)-P,0) 

87215333 

33,3. 



H(*,A,MA) c T3(1) 

87215334 

33«. 



H(NA,EU) « T3(2) 

87215535 

335. 


730 

H(EN,NM ■ 0.0 

87215336 

336. 



H(EA,FN) « 1.0 

87215357 

317. 



ENM? 9 NA . t 

87215338 

33fi. 



IF.(EN'^2 .E(J. 0) GO TO "on 

87215359 

339. 

c 



87215340 

3«0. 



no 790 11 9 1, ENP? 

87215341 

3U1. 



I 9 NA - TI 

87215342 

342. 



*> = H(I,I) - P 

87215343 

343. 



PA s 0.0 

87215344 

344. 



SA * H(I,EN) 

87215345 

345. 

c 



87215546 

346. 



DO 760 J 9 M, Mfl 

87215347 

347. 



RA 9 PA ♦ * H(J,NA) 

87215348 

348. 



5A 9 SA ♦ • H(J,FN) 

S7215349 

349. 


760 

rONTINIIf 

87215350 

350. 

c 



87215551 

351 . 



IF (wj(l) .GE. 0.0) GO TO 770 

87215552 

352. 



71 = I'l 

87215353 

353. 



P s PA 

87215354 

354. 



S s SA 

87215355 

355. 



GO (0 790 

87215356 

556. 


770 

P 9 T 

87215357 

357. 



IF (Wim ,NEc n.O) GO TO 7"0 

87215358 

358. 



73 = CMF'LX(-RA,-SA) / CMPLX(W,0) 

87215359 

359. 



H(I,MA) 9 T5(l) 

87215360 

360. 



H(1,EN) s T3(2) 

87215361 

3ft 1 . 



, GO TO 790 

87215362 

362. 

c 



87215363 

363. 


780 

X = H(I,Ifl) 

87215364 

364. 



V » H(I+1,I) 

87215565 




39 


385. 



V8 ■ (INRC!) - 8) • (W8(n - P) ♦ Wld) $ 

NKI) -0*1 

388. 



VI • (NR(n - P) 4 2,0 # 0 


387. 



IF (VP ,E(9. 0.0 .AND. VI .CO. 0.0) VP ■ 

MACHCP 8 NQPN 

368, 



* • (A88(W) ♦ A8S(0) ♦ A8S(V) ♦ A8S(V) ♦ 

ABS(ZZ)) 

389. 



23 « rHP(.X(V*P-ZZ*PA«QftSA.X*S«2Z4SA-08PA) / CPPLXtVi.' 

370. 



H(I,NA) ■ T3(l) 


371. 



H(|,£N) m T5(2) 


372. 



IF (ABS(X) ,LC, A8S(ZZ) ♦ A8S(0)) GO TO 

785 

575, 



H(I41,NA) ■ (*04 • «< * H(I«NA) 4 0 • H(I 

♦CN)) / X 

374. 



H(I41,EN) • (-SA - Kt • H(I.CN) - Q * H(I 

»NA)) / X 

575. 



CO TO 790 


578. 


785 

Z5 « CMPLX(-R-y*P(I,NA),-S-Y«H(I,EN)) / 

CPPLX(ZZ,0) 

377. 



H(U1,NA) ■ 73(1) 


378. 



H(I*1,FN) » T3(2) 


379. 


790 

CONTINUF 


380. 

C 




i«l. 


80S 

CONTINUF 


382. 

C 




583. 

C 



* 

384. 



no 840 I * 1, N 


385. 



IF (I .OF. LOW .»Nn, I .Lt. IGH) GO TO 840 


388. 

r 




387. 



00 821) J « I, N 


3«8. 


820 

Z(1,J) * 


389. 

c 




390. 


840 

CONTINUE 


391. 

r. 


A********* multiply by TWANSFORMATION »'ATRIX TP GIVE 

392. 

c 


VFrTPPR PF PPICINAL FULL P8TR1X. 


303. 

c 


FOR JtN STEP -1 UNTIL LP8 DO -- *** 


3914. 



Pp 880 JJ 1 LPHi, N 


595. 



J e Ni ♦ LOW • JJ 


598. 



M t MjNOfJ.K.H) 


397. 

r 




39fl. 



DO 880 I * LOW, IGH 


399. 



ZZ » 0.0 


400. 

c 




401. 



no 880 K S LOW, M 


402. 


880 

ZZ « 7Z ♦ Z(I,K) • H(K,J) 


405. 

c 




40(1. 



Z(1,J1 * 7Z 


4 05. 


880 

cnuTtNiif 


408. 

c 




407. 



RO T() 1001 


40H. 

r 


*****AA«#A SFT fpRPR -- NO CPNVFHGFNCF TR AM 


409. 

c 


fclGFNVALUf after 30 ITFRATIONS **** 


410. 


1 noo 

IFRR s FN 


411. 


tool 

RETURN 


412. 

c 




41s. 



END 



S72I5S** 
BTZiSSfcT 
A72I9SM 
B72l5S6f 
B72I5S70 
»721«7I 
B7215S72 
S72JSS7S 
B7215S7il 
B7215S7*5 
S7215S7B 
B7215177 
«721'»S7(I 
872t5379 
B7215JB0 
S72l5S8t 
872I5I82 
87215J83 
87215584 
87215385 
87215388 
87215387 
87215388 
87215389 
87215390 
87215591 
87215592 
87215393 
87215394 
87215395 
87215398 
87215597 
8'’?15398 
87215399 
87215400 
87215401 
87215402 
87215405 
87215404 
87215405 
87215408 
87215407 
87215408 
87215409 
87215410 
8721541 1 
87215412 
87215413 
87215414 
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{ 

1 


1 . 

2. 

5. 

a. 

$. 

6 . 

7, 

10 . 

II. 

ta. 

13. 

10 . 

15. 

Hi. 

17. 

1^. 

10 . 

20 . 

? 1 . 

22 . 

23. 

20 . 

25. 

2 <). 

27. 

2H. 

20 . 

30. 

31. 

32. 

33. 
SO. 

35. 

36. 
57. 
3«. 
50. 
00 . 
01 . 
02 . 
03. 
00 . 

05. 

06. 

07. 

08. 
00 . 

50. 

51. 

52. 


SUBROUTINE BALROK (NH,N.LOM* ICH, SCALE «M, Z) 

C 

INTEGER 

REAL 8CALr(N)»Z(NN,M) 

REAL S 

C 

C THIS SliflRnUTlNE IS A TRANSLATION Of THE ALGOL RROCEOUiE BALBAB# 

C NIJM, math, 13, 203-300(1969) BY RARLETT AND REIN8CH, 

C HANOBOOM FOR AUTO, COMP,, VOL . I 1-L INCAR ALGEBRA, 315-326(197} ) . 

C 

C THIS SUBROUTINE FORMS THF FIGENVECTORS OF A REAL GENERAL 
C MATRIX BY BACK TRANSFORMING THOSE OF THF CORRESPONDING 
C HALANCFP Matrix OETERMINEO by balanc. 

c 

C ON INPUT- 

c 

C NM MOST HE RET TO THE Rnv>' OIMENSIOM OF THO-OIMENSIONAL 

C ARRAY PARAMFTERS AS I'FTLAREO IN THE CALLING PROGRAM 

C dimension statement, 

C 

c N IS the order of the matrix, 

c 

r LOW AND IGH APE INTFGERS DETERMINED BY RALANC, 

C 

C SCALE rONTAINS TNFORMAT ION DETERMINING THE PERHUTATIONS 

C ANO sealing factors USED BY BALANC, 

C 

C M IS the NyMBER OF COLUMNS OF I TO BE BACK TRANSFORMEO, 

C 

C Z CONTAINS THF REAL AND IMAGINARY PARTS OF THE EIGEN- 

r. VETTnRS TO BE HAfK TRANSFORMEO IN ITS FIRST m COLUMNS. 

r 

C ON fliiTPUT- 

C 

C 2 CONTAINS THF REAL AND IMAGINARY PARTS OF THE 

C TRANSFnRMfD EIGENVECTORS IN ITS FIRST M COLUMNS. 

C 

r OUFSTIONS and comments SHOIILO be OIPECTEO to B, 8. GARBOW, 

C APPLIED maTHFMATTCS division, ARGONMf NATIONAL LABORATORY 
C 

C ......... — 

C 

IF (TGH .FD. LOW) GO TO 120 
C 

on IK) I 3 LOW, TGH 
S 3 SCALEd) 

C ,«*****«** lEFT HAND EIGENVECTORS APE BACK TR«NSF0HME0 

C IF THF FOREGOING STATEMENT IS REPLACEO BY 

c s»i .0/SCALE ( n . **#•****** 

’ DO ion u s I, M 
KiO Z(I,J) * Z(I,J) • S 
C 


70215009 

70215005 

70215006 

70215007 

70215008 

70215009 

70215010 

70215011 

70215012 

70215013 

70215019 

70215015 

70215016 

70215017 

70215018 

70215019 

70215020 

70215021 

70215022 

70215023 

70215029 

70215025 

70215026 

70215027 

7021502B 

70215029 

70215030 

70215031 

70215032 

70215035 

70215039 

70215035 

70215036 

70215037 

70215038 

70215039 

70215090 

70215091 

T0215092 

70215093 

.70215099 

70215095 

70215096 

70215097 

70215098 

70215099 

70215050 

70215051 

70215052 

70215055 

70215G59 

70215055 




SI. 


no 

CONTINUe 

70215056 

54. 

C 


**•••***. roR IBI.ON-] STER -i UNTIL 1# 

70215057 

55. 

C 



70215050 

56. 


120 

DO 140 IT ■ 1, Ni 

70215059 

57. 



I ■ I! 

70215060 

5«. 



IF (I .Gf. LOW .AND. I .LE. TGH) CO TO |40 

70215061 

59. 



IF (I .IT. LOW) I « LOW - II 

70215062 

60. 



K ■ SCALE(I) 

70215063 

61. 



IF fK ,ro. I) GO TO 140 

70215064 

62. 

C 



70215065 

63. 



DO 130 J * 1, M 

70215066 

64. 



5 « Z(I,J) 

70215067 

65. 



7(1, J) » Z(K,J) 

70215068 

66. 



Z(K,J) s S 

70215069 

67. 


130 

CONTTNUF 

70215070 

66. 

C 



70215071 

69. 


140 

CtiNTINUfc 

70215072 

70. 

C 



70215073 

71. 



f^ETHKM 

70215074 

72. 

c 



70215075 

73. 



I^ND 

70215076 
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1. SUMOUTINE CISORT(WB»Wl,rH,TC,ZETA»OMC,Pl,N) 

2. 01*’EN»ION WR(n»Win).FP(n»TC(n,2FTICJ),SB(5«)»SI(50) 

5. C 

C s« *NO SI »PF TEMPOPARV WORKING STORAGE 

5. C 

6. JUMRil 

7. NT»N 

8. MRBO 

R. NV1*N-t 


10 . 


no 100 i«l,u 

n. 


lF{JtlMP)?n,tO,2C 

12 . 

10 

JIJMP«1 

n. 


GO TO 100 

la. 

20 

!F(Wl(nH0,10,60 

IS. 

30 

NR*NR4l 

18. 


SR«RR(I) 

17. 


no SO 1 * 1 .NM 

ta. 


IF(L-n 50 , 90,/|0 

IR. 

90 

WR(L)«k*R(lH) 

?0. 


i^KI )«1'.I(LM) 

21 . 

SO 

CONTTNl'F 

22 . 


WR(NT)rS'^ 

23 . 


>«-I (NT)«0. 

2 «. 


NTsfvT-1 

2 *^. 


f,n TO too 

2h . 

60 

,r.|MP*0 

27. 

ion 

CONTINI'f 

2 H. 


NMPrN-^R 

2 R. 


NMR1*N^'R-1 

50. 


00 200 I« 1 ,NMR 1 ,2 

51. 


W1I«WI(1) 

52 . 


*»TR»wP( t ) 

55. 


IWAKsI 

3«. 


lP2aIt? 

3 S. 


1F(IP2-NMR)1U0, 160, 160 

38. 

190 

on ISO J«IP 2 ,n»'R ,2 

37. 


IF(wii.(;r..6i(j))r.o in iso 

38. 


wn«wi(j) 

39. 


'^lRi>*<Rf J) 

90. 


IMAX=J 

ai. 

ISO 

CONTIMIIF 

« 2 . 

160 

SRdlawRU^'AX) 

U3. 


SR(I*n=rtR(lt'ftX) 

un. 


SI (I)e'Al (TMAX) 

«S. 


SI (!♦! )*-wi UMAX) 

96. 


WR(!M4X)swR(I) 

97. 


wH(IM4K+n=W9(l) 

9H. 


WI (IMAXlsWI Cl ) 

99. 

■ 

WI (IMAX^l )«-Wl Cl ) 

SO. 


ZtTACl )=SR(n/SQRT(SRn)*SR(I)FSI (n*SI(in 

SI. 


ZFTACltl )=ZETA(I) 

52. 


TC(I)«l./(2.*Pl»SR(n) 


53. 


TCn»i)«TC(!) 

5«. 


r«fn«»i(n/oMG 

55. 


fR(ltU»FH(n 

56. 

299 

COMTI^UC 

57. 


IFf5R)2t9,260,?l9 

55. 

210 

90 250 r«NRB,N 

59. 


sR(n«w9(n 

69. 


91(1)«61(I) 

61. 


7ET6n)»l. 

63. 


Tcn)*i./(?.*pi*5H(m 

63. 


FR(n»o. 

6U. 

250 

CONTINUE 

65. 

260 

DO 500 X«1»N 

66 . 


WR(l)«8F(n 

67. 


NI(l)*5I(n 

69. 

300 

CONTINUE 

69. 


RETURN 

70. 


ENO 


APPENDIX B 


DEFINITION OF SYMBOLS 
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.SYMBOLS 


ABO 

ADUM 

AKA 

ALA 

AKAO 

ALAO 

AKB 

ALB 

AKBO 

ALBO 

ALPHA 

ALT 

AM 

APHI 

AREF 

B 

BKA 

BLA 

BKB 

BLB 

CLA 

DELPSI 

DEN 

DI 


Aero term 
Dummy array 
transformed aero term 
Aero term 

Transformed aero term 
Aero term 

Transformed aero term 
Aero term 

Transformed aero term 
Aero ferm 
Angle of attack 
Altitude 

Aerodynamic torque 
Initial transformation mati'ix 
Refei'cnce area 
Matrix 

Transformed aero term 
Aero term 

Transformed aero term 

Aero term 

Aero term 

Delta PSI 

Aero term 

Dummy index 
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SYMBOLS (Continued) 


DREF 

BUM 

EUM 

FR 

FV 

GG 

GM 

lERR 

IV 

LI 

LIMIT 

MI 

Ml, M2, 

OMG 

PAMA 

PAMP 

PAMS 

PAS 

PASI • 

PGGS 

PHI 

PI 

PHIB 


Reference diameter 
Dummy array 
Dummy matrix 
Frequency of root 
Working array 

2 

Gravity constant, 3H 
Universal gravitational constant 
Ei’ror flag 
Working array 
Woi’king storage 
Number of iterations 
Working storage 
M3 Aero terms 

Orbital angular velocity 

Partial of aerodynamic torque with respect to alpha 
Partial of aerodynamic torque with respect to PHI 
Partial of aerodynamic torque with respect to PSl 
Partial of alpha with respect to PSI 

Partial of aerodynamic torque with respect to PSI divided 
by momeht of inertia 

Partial of QX li’ ~ gravity gradient torque, with respect to 
PSI expressed in the principal axes reference frame 

Aerodynamic roll angle 

tr 

Aerodynamic bank angle 
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SYMBOLS (Continued) 


PNUS 

PNUSI 

PPHI 

PPS 

PSI 

PS 10 

QI 

QNUS 

RAD 

RE 

RHO 

RHOK 

RORB 

TC 

TCP 

TMM 

TPG 

U 

V 

VMAG 

VO 

WI 

WR 


Partial oC NU with respect to PSI 

Inverse of the partial of NU with respect to PSI 

Initial isolation angles 

Partial of PHI with respect to PSI 

Rotation angles 

Initial input of rotation angles 

Inverse of angular velocity transformation matrix 

Dummy matrix for PNUS 

Conversion factor radians to degrees 

Radius of the Earth 

Atmospheric density (slugs per cubic foot) 

Atmospheric density (kg per cubic meter) 

Orbital radius 

Unstable time constant or logaritmic increment 
Transformation matrix from principal to geometric system 
Characteristic polynomial matrix 

Transformation matrix from geometric to principal system 

Transfoi’mation array for converting aerodynamic torques 
from geometric system to principal system 

Velocity in geometric system 

Magnitude of velocity vector 

Velocity in principal system 

Array containing imaginary parts of roots 

Array containing real parts of roots 
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SYMBOLS (Concluded) 


XCG X coordinate of center of gravity 

XI Moment of inertias 

XNU NU 

XNUMAG Magnitude of NU 

YCG Y coordinate of center of gravity 

ZCG Z coordinate of center of gravity 

ZETA Ratio of the real part to the x^oot sum square of the real 

and imaginary parts of the root 


49 


APPENDIX C 


INPUT DATA 


TABLE 1. INPUT REQUIREMENTS 



Ib’OHJ’niii 

Hyiidjolh 

1 

DU HI- 
AHUI* 

•> 

M 

UMIT 

;} 0 

TiilHl..!) 
1 i.:i 
,i i.;t 

li 

Xll 


XI2 


xi;i 

7 

NC’ti 
YC’U 
. zed 

S 

Ml 

M2 

m:i 

i) 

A 110(1) 
1 l.Ml 

10 21 

AbAOd..)) 
Al, 110(1 ..1) 
.) 1.M2 
1 1 .M 1 

22 1G5 

Ai,A(l ..I .K) 
llbAd.o.K) 
K 1 di:i 
.1 l..vi2 

• 

Al.lld -I .K) 
111,11(1 ,.1 .1\) 
K i.m:i 

,1 l,Ai2 


1 l.Ml 

KHi 

HIlOK 

AI,T 

lf)7 

HHIOd) 


I l.:i 
ir.K 


Vni'iiibU* DcrinUlon 

UflVi'tMiL'o 1,1‘nn'th 

Ut'lOlH'IU'O Al'C'II 

Numhoi* of lUM'iilioiiH 

Ti'anKrurnmUon MiilHx 

Iroiii Pi'liu'i|)i»l lo 

(U’oiJU’li'ii* C’tMn'(iIiii)U'.s 

X ainincnl of 

huM’Un 

V I'lMlifipiil Moiiiom «»r 
liu'Hiii 

Z I'l'lui’ipnl .Mtmu'lU ol' 
hu’i'dii 

(U't)iiu’l rii' C'nortrmaU's 

ol iX'iUoi* oT (li'iivity 


Aoi’otlyiuiink* Tfiblo 
liulicics 


(’ooirUMonth or 
l)otil)lc’ l*(nu'k'i’ SoiMOh 
Ibi- Ai'i’odyiiiimif FoI'cph 
timl Toi’rivic'H 


AtimisphcM'io DoiisUy 
Allilutk* 

An^rU' of Allm’k . 

Hull AniJla. Hniik 
Anp;lo 

Hl.-iiik Cm'd 


UnilK Formitt 


Ft. 

0 

K(|. Ft. 

(i 

IJnltlCKH 

I 

IblltlOSK 

30 


Klu{? St|. 

Ft. 

G 

Hluu Ut|, 

Ft. 

G 

HluR S(i, 

Ft. 

G 

Ft. 


G 

Ft. 


G 

Ft. 


(3 

Unitloss 


Ib 

UnitlpHH 


15 

tInitlOHH 


15 

dnitloKH 


F12.0 

UnitlosH 


F12.0 


Kp;/cu. 111 . 

G 

Niiutioid Milo 

G 

Dogroes 

3G 
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?0 

.l70O^Om,-,2')7^'nifc»,.9( 
PI. t'lawi.ri. . ^aWVuOKCY. 


H - 

6 6 6 


9 - 

.1)00 

,nno 

in - 

7.8a? 

,npn 

n - 

-.125 

.non 

i? - 

-.201 

.OOP 

1^ - 

.059 

, (1(1 II 

H - 

.1)15 

,nr.p 

1 *> - 

-.il6’l 

,0(111 

U> - 

-1 .?n9 

,oon 

17 - 

-.057 

.non 

in - 


.(•00 

19 - 

.166 

. 0 ('n 

2n - 

-,''56 

.non 

21 - 

. »51 

.nfi’) 

2? - 

.nno 

.727 

2^ - 

.ni'* 

• . 1 2n 

?« - 

.''111' 

.non 

25 • 

.■nc 

. '! n n 

2h - 

'06 

• # 

-.532 

27 - 


.115 

2H - 

. inn 

.non 

29 - 

.mill 

.fun 

in - 

,'UJO 

.nnn 

ii - 

,l'Ou 

.ilflO 

12 - 

.'I'lfi 

.np( 

11 - 

.iiiiq 

.non 

in « 

- , 260 

-.010 

55 - 

-.'111 

.021 

16 - 

.non 

.onn 

17 - 

,000 

.no( 

IH - 

-. 2 'in 

-.niiR 

59 - 

.001 

.on 

(in - 

.00(1 

.nno 

'11 - 

.i>nn 

.non 

«2 - 

.i'OO 

.onn 

ui - 

. "00 

.OP'i 

«i/i - 

,"on 

.non 

'in - 

.linn 

.imn 

(i() - 

-.07? 

-.171 

UI - 

.’11 

« , n?2 

HR - 

.20I“ 

. 0 n 11 

no - 

,'UIO 

, (Hill 

sn - 

.■10 7 

.109 


-,|7t 

.?87 

,A?ii 

-.127 

-,nj9 

,i*n? 

-,nu/l 

,<*7P 

.('2b 

.f*n 

,f At 

.i'?h 

.I'nn 

.llAA 

,aa.. 

.tin'' 

.on- 

-.o?f 
.()('(' 
.no I: 

-.o\7 

-.IIJ6 

.nno 

.Mill 

.llO( 

.(ift'- 

.nftt 

,iin/ 


HfcFfRtNCe 0|*Mf7E» ^HO I»< rT»FT2 

NUMBFB of ITE»*T!0N8 
TFAHSFOFM/kTION FROM 
FFIHflPAL TO CEOMETMIC 
COnROlNFTFS 

PRIMCIPIL INERTIAS IN SLUO-FT2 
(OfATION OF CO IN FT 


non 

,000 

.OOP 

,000 

.000 

.000 

,000 

• OOP 

.000 

,000 

,000 

.00(1 

.11 no 

.000 

,000 

,000 

,pon 

.000 

, 0 On 

,000 

.000 

,(I0(' 

.OOP 

,000 

,000 

,000 

,000 

,0P() 

,000 

.000 

,.100 

,00(1 

.000 

.('OO 

,()np 

.OltO 

.oon 

.000 

.000 

.076 

.01 a 

.107 

."3? 

-.012 

- , 006 

,i;oo 

,ii()(' 

.OtlO 

,1)00 

.1)00 

.non 

.071 

-.002 

-.021 

.032 

,019 

.005 

,000 

.OOP 

*000 

.000 

.non 

.000 

.non 

.OOP 

.<100 

.DPP 

.000 

.000 

,npp 

.OOP 

,000 

.100 

.000 

.000 

.OOR 

,013 

-,P16 

.00 7 

-.005 

,0P1 

.000 

,000 

,000 

,( on 

,000 

.000 

.015 

.nil 

-.016 

,007 

.00(1 

,000 

,1)0(1 

.000 

.000 

.non 

,000 

•lOOP 

,pon 

.onn 

.000 

.000 

.OOP 

.OOP 

,000 

,000 

,000 

,000 

.non 

,000 

-. )?3 

-.007 

.0 30 

• , 0 1 fi 

.005 

-.010 

.non 

.onn 

.000 

.ilOO 

.poo 

.opo 

.03H 

,011 

-.039 
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51 - 

• , U 1 6 

.036 

-.020 

,027 

-.005 

.016 

52 - 

. 1)30 

.000 

.000 

• OPO 

.000 

.coo 

■33 - 

. I'OO 

, 90 u 

. 1/00 

,000 

,000 

.000 

5 « - 

.990 

.000 

.000 

.000 

,000 

.OPO 

55 - 

.090 

.009 

.900 

.000 

.OOP 

,000 

5 A • 

. 1)00 

.000 

.000 

.000 

,ppp 

.000 

57 - 

.000 

.009 

.000 

.000 

,000 

.000 

5 f > - 

7.296 

.005 

.679 

-.016 

-, PP 2 

,005 

59 - 

-.162 

.027 

i *, 10 ? 

-.002 

.009 

-.007 

f )0 • 

,000 

.000 

,009 

.000 

,000 

.000 

ftl • 

.990 

.000 

, l'OP 

,000 

,000 

,000 

- 

- 1.079 

. 00 ) 

-.156 

-.OOP 

.023 

,000 

63 - 

. Jl ? 

.015 

.011 

,003 

-.022 

,000 

bit - 

.990 

. i )0 9 

.OOP 

,000 

.000 

.000 

65 • 

..jOn 

,900 

.909 

,000 

.poo 

.000 

66 “ 

.100 

. 90(1 

,990 

.OPO 

,() 0 P 

.OOP 

b 7 " 

.000 

.OOO 

.(•on 

.000 

.000 

.000 

6 « - 

,909 

, 9 nfi 

.OOP 

.000 

.000 

.000 

69 ■ 

.900 

.(■no 

.'lO'i 

.000 

,000 

.000 

76 - 

-. 1)1 

.906 

,007 

-. 01 ? 

,003 

. 00 « 

71 - 

.009 

.297 

.011 

.039 

’ -.011 

.017 

7 P - 

.(•90 

.000 

.009 

,noo 

.000 

.000 

71 - 

. 00(1 

, 0 0 1 ; 

,900 

,000 

.non 

.000 

70 - 

,' M 6 

-.996 

-.012 

,015 

-.002 

-, 005 

75 - 

-.0 1.7 

-.099 

- . 0 1 ,1 

-.955 

.000 

-, 0 ? 5 

76 - 

.0 9 .) 

.non 

.non 

,000 

.900 

.000 

77 - 

' .910 

. 00(1 

,9011 

,000 

.OOP 

.000 

76 - 

.1(011 

.0911 

. 0(10 

.000 

.000 

.900 

70 . 

.|| 0 ([ 

.(■or; 

, 1)09 

,1100 

,000 

.009 

rtn • 

. o ‘)0 

,090 

.009 

,000 

.090 

.009 

Ht - 

. XKl 

,900 

. (1 0 (1 

.001) 

.009 

.000 

- 

.) S 0 

.025 

-.206 

.027 

.050 

-.031 

«3 - 

-5.274 

- , (196 

-.925 

,()«a 

.196 

-,000 

HO > 

.900 

.900 

,000 

,000 

,000 

.000 

H5 - 

,0011 

,900 

,nno 

.009 

.OPO 

.000 

66 - 

-.029 

.007 

.059 

.non 

-.019 

.001 

67 - 

.7)6 

-.014 

.225 

.non 

-.050 

.000 

K6 . 

.009 

. (’OO 

.999 

.900 

.000 

,000 

HO • 

.090 

,900 

,( i 09 

.000 

.000 

.000 

90 - 

.90.) 

. 00(1 

,900 

.900 

,900 

.900 

91 - 

. n.i 

,(,00 

.009 

.000 

.nOo 

.009 

9? - 

.1)00 

.non 

.090 

,000 

.non 

,000 

93 - 

.((Of. 

.000 

.0 09 

.000 

.900 

.000 

90 - 

.1 55 

-.194 

-,nno 

.080 

.007 

-.037 

OS - 

-...O') 

.009 

-.01 ) 

.('22 

-.005 

.011 

96 - 

.OOll 

.090 

.non 

,000 

.POO 

.000 

97 - 

. ('0 (1 

, 9 9 (' 

.non 

.000 

.000 

.OPO 

96 - 

-.1173 

. 060 

.(inn 

“ , 0 7 6 

-.011 

.033 

00 • 

, 9 n 9 

-.9?2 

,014 

-.036 

.007 

-.020 

1 OO - 

.000 

. 0 0 0 

.0 9 0 

.900 

.900 

.000 

16 ! - 

.009 

.090 

.99 9 

.OOP 

,000 

.090 

1 '1? - 

.non 

.non 

.0 09 

.000 

,000 

.OOP 


ORIGINAL PAGB is 
ROOR QVALI'JY 
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lOS - 
loa - 
105 - 
t06 - 
107 - 

I Ofl • 

|()o - 
no - 
n I • 
n? - 

iij - 
lio- 
ns - 
n 0 » 

117 - 

II H . 

no - 
121 - 
i?i - 
12? - 
l?J - 
l?o - 
1?S - 
l?6 - 
127 - 
I?** - 
129 - 
ixn . 
Ill - 
13 ? - 
135 - 

no - 

I3S - 
I 36 « 
157 - 
I3« - 
I 50 - 
I '10 - 
Mil - 
Ml? - 
Ml 5 - 
I '10 - 
lOS - 
1 06 - 
107 - 
Ml« - 
100 - 
ISO - 
1S| - 

n? - 

|S5 - 
ISO - 


,000 

.000 

,Ono 

,000 

.nno 

,000 

-5,106 

.01! 

.100 

-.017 

. non 

.OflO 

.000 

,000 

.779 

-.007 

-,P 10 

-.01? 

,i)00 

.000 

,"00 

,000 

. vll/(t 

.000 

, lOf, 

.000 

, non 

,00( 

,0"tl 

.non 

-,"67 

-. ?66 

-,O0h 

-.OSS 

, non 

,0"0 

, iir.i. 

, "(.<1 

, 116 

.?0 7 

.OSS 

.031 

, I'Hl) 

.00" 

, MOn 

,uon 

.•Ifi'l 

,00" 

.non 

.orn 

,"0 0 

.0(0 

,00"' 

.non 

- , 1 0 i| 

-."Ml 

J . ? 1 n 

.0?? 

.no 1 

.00/) 

,ono 

.{00 

.07? 

.003 

-.067 

. 0 0 '1 

.O0(i 

.00" 

,n n 

. iMiO 

,00 0 

, "Of, 

. jO'I 

, 0 n 0 

,00" 

,oro 

. '"0 

. non 

-,ii7t 

-.nJH 

, i‘. 1 rt 

-.059 


.00" 

,t 0‘l 

.u"'l 

.007 

, nsn 


.'in? 

.00,, 

.'100 


.00,1 

. 0 n n 

.00" 

,000 

.non 

. f 0 " 

. ono 

. 0 '■ a 

.0",) 

' . 1 3 n 

.0?,| 


.900 

,000 

.|■’l{> 

.000 

,"00 

.900 

-.350 

.OMl 

,l?ii 

.00? 

.00(. 

.000 

.000 

.(>00 

,l?3 

.005 

-.')?l 

.001 

,000 

.000 

.noo 

.000 

.000 

.000 

, 1 0 0 

,000 

, 0 0 .1 

,000 

.000 

.00 0 

.n|o 

.016 

n?6 

-, oso 

.no,, 

.non 

.1.00 

."00 

-. 0"7 

-.03? 

. n Ml 

.068 

, "00 

.nof) 

."00 

,000 

.1100 

.000 

.00,1 

.000 

,000 

.non 

."Oil 

.000 

. ?60 

-.00? 

.ns'i 

-.OS? 

.rno 

,000 

.0,)0 

.000 

-.360 

-.115 

-.130 

.11 0 

.0 no 

.,100 

. >‘0" 

,000 

."0" 

.000 

.'ill 

.000 

.00" 

.non 

,n 0" 

,00f 

. (iiiO 

-."1? 

.'M3 

-.117 

,0(1 0 

.100 

.r.0( 

."00 

-.""9 

.11" 

- . '1 1 6 

.o?8 

."00 

.000 

.00 1 

.10 0 

,0 00 

.100 

.00(. 

.000 

,000 

.000 

. n (1 ('■ 

, non 

.?S0 

-."OS 


.000 

.000 

.000 

.000 

,000 

.000 

-.010 

-.007 

-.065 

.000 

,000 

.000 

,000 

.000 

-.03? 

-.001 

.03? 

-.001 

.000 

.000 

.000 

.000 

.00(1 

.0(10 

.noo 

.000 

,000 

.000 

.000 

.000 

-.006 

,01? 

.018 

-.009 

.000 

.000 

.no(' 

.000 

.006 

-.OH 

-,0?7 

,01? 

.00 0 

.000 

.non 

.011' 

.000 

,000 

.000 

• oon 

.non 

.000 

.000 

.000 

-.070 

.006 

-.055 

.003 

.000 

.000 

.000 

.000 

.05? 

.000 

.003 

-.001 

.000 

.000 

,0"o 

.000 

.non 

.000 

.000 

.00(1 

.000 

.000 

.non 

.000 

-.101 

.006 

.00? 

-.011 

.non 

.000 

.1')C 

.001) 

.nos 

-.016 

.noo 

.018 

. nnr 

.000 

.non 

.one 

.ono 

.000 

.000 

.000 

.000 

.000 

.000 

. noo 

-.050 

-.007 
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iss - 

,6U1 -.016 

isfe " 

,000 ,nc<i 

IS? - 

,000 .00(1 

l'5H - 

-.210 .006 

tS9 - 

-.162 ,1105 

1 ” 

,00(1 .00(1 

- 

,l)0(^ .OOA 

162 - 

,:no .non 

16^ - 

,000 ,000 

16/1 - 

.non , 0 (,o 

16S • 

, ( 1 0 0 . 0 n n 

1 66 ■ 

.?26«;f-11,?55. 

|67 - 

1S9.??1 ,81 .266,0.666 

166 - 

, 266 E,-n , 250 . 

166 - 

156<,6io,pl .501 ,5.0/iS 

170 - 

.1065F-1 1»??5. 

171 - 

150.618,81 .755,5,16(1 

172 - 

.161F-1 1 ,220. 

17< - 

15?.m-27,H2.0?n7,5.72185 

17'4 - 

.0P15F-1 1 ,?15. 

17S - 

ll7.,rt/l.,o. 

176 - 

.06Jt-n ,210. 

177 - 

1 y8. 17A,H2.5b7,6.30 S 

17« - 

.A7A5P-1 1 ,?05. 

176 - 

1 OO . , AO .,7,6 

1 >\n - 

.6P'l5f-n ,200. 

Irtt - 

Ml ./l,ft5.P#7.?o 

IP? - 

.flOlOt-11 ,165. 

IPl - 

l'i(l.,AO., 7.5 

IH'l - 

.607sr-1 1 , 1°'i. 

105 - 

116., 60., 7. 7 

1P6 - 

.n?85F-l/',lA5. 

1H7 - 

10 'J, 6^2, 81.6^1, 7, ?61 

) AA . 

.tlin5p-i(i,lA0. 

1 AO - 

1 IP. 2, "5. 1!, 7. •VI 

160 - 

.16nnF-lO,175. 

161 - 

1 16.0,85. ,7.7 

)0? - 

.l«25F-10,17n. 

161 - 

138. 5, 85., 7. 7 

1 60 - 

.2115F-1 1,165. 

16>; - 

1 17.626,65.116,7.71 \ 

16b - 

.276SF-!0, 160. 

167 - 

tl6.?,rt-j.n,7.6 

lOA - 

.l(|„r,F,to,155. 

166 - 

1 lb., 87., 7. > 

200 - 

.Olhiil -10, 15 1. 

?l'l - 

116. ,88. ,7. 6 

?n? - 

.51 2K-1 1,105. 

?C1 - 

1 IS, 7n1 ,87.515,7.680 

- 

’.6l0f-1i>,100. 

?06 - 

135.5,80.(1,8. 1 

206 - 

.761SF -1 0, 1 1S. 


050 

,005 

000 

.000 

000 

.000 

055 

-.005 

096 

,000 

000 

.000 

no(‘ 

.000 

non 

.000 

ron 

.000 

got' 

.000 

onn 

.non 


052 

-,005 

000 

,000 

000 

.000 

019 

-.003 

007 

-.002 

OOO 

.000 

000 

.000 

000 

.000 

lOOO 

,000 

.000 

.000 

,000 

.nco 
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207 - I JO.«,90,S,O.2 

209 - ,99a«;r-io,n«. 

2f)«> - 15". 6, ot .ft, 0,9 

210 - .12t>(it-0,|?s, 

211 - 130. 9,92, 9, ft, 5 

212 - .Ife2fc-9,i20. 

213 - 13". 2, 99.1, ft. 7 
219 - .?10F-9,n5, 

215 - 139,127,99.57 J, ft. ft9»i7 
21h - .?7^>^-9,| 10. 

?17 - 133. ft, 97. 5, 9. 5 
21ft - .36HftP-9,105. 

219 - 133.7fto?9,9fl.i«95,9.7599 

220 - .S0l5f-9,ton. 

221 - 133.ft3(j,inn. 197, i,).53('3 

222 - . 6991 - 0 , 95 . 

223 - 133. 519, 101. 975, 11. 65ft 
229 - .tnO05F-ft,9O. 

225 - 133.0, 109., 1 l.ft 

226 - . I5('ftb--fl,ft5. 

??7 - 153.353,105.353,16.191 
22fl - .20!ip-fl,90. 

229 - 133. 30P, 106. 7&5, 21 ,17ft 

230 - 

231 - 

232 - 




APPENDIX D 

PROGRAM OUTPUT AND PLOTS 


ALTITUOe 21S 

PHr 

.l2l«30f4E-l6 - 
HEAL 

-.IHfc90156E-OA 
-,lP690|56E-nA 
920259 lSe-(, 6 
-.920259ISE-O6 
,#in713570E-0(i, 
,^071 3570f-'tA 


ALTITIIDE ?30 
PHT »59.9IO«7h 
.l21'450t)aE-l^ 

PE»L 

•.579SS«e?F«06 
-,3795Sfi«2E-06 
-.5?693?71E-''h 
-.5?iS9327lE-( h 
.9('«59<J1S3F-o6 
.OOS'»9I*53E-06 


ALTITliDE 22‘> 
PHI IS9.M7S)9 
-.23«lH767E-t6 

heal 

-,5nSt70fl3E-''6 

-,5PSl70fl3F-''^ 

-.6Sfi6P9<>2E-nfc 

-.6S069962E-''iS 

,1?3'5«65^F-0S 

.laSbH^SS^-ns 


ALTlTtine ??n, 
PHI lS2,9n?6<>-5 
-.19')9P32nf-if, 

RFAL 

-.7A755729F-16 
-.7H7‘557?«F-.)ft 
-.70900 |93F-f.A 
-.799fl0103F-r6 
. IS«19792F-f'5 
.tS'n0792F- IS 


.0 OEHSITV .229S0E-H 
alpha P1.269A9A PHIfl 9.6S6SP9 
.S77662V2E-1S -.5209 1 709f -I 7 .S77P1397E-I5 


JMAO 

FHFO 

2ETA 

TC 


.2079A799F-02 

I.A527S 

-.0001 

-.A5E 

06 

-,?079«790E-02 

I.P5273 

-.0001 

-.«5E 

06 

.17915523F-0? 

1.55135 

-.0002 

-.5<IE 

06 

-.17/M5525F-02 

1.55135 

-.0002 

-.38E 

06 

.9319A157E-03 

„3A939 

.0019 

.26E 

06 

-,9n96I57E-03 

.30939 

.0019 

.26E 

06 


0 I'FNSITV .?660(lf-n 

alpha Pl,SCI26f< PHIH S, 095900 

9391 'IA59F-I*: - .2? UlO 377F- 1 0 .23993899E-I 9 


THAR 

FHFO 

7FTA 

TC 


.?')*136125E-0? 

t.85??7 

«.0002 

-.9?E 

06 

-.?0H36l?SF-O? 

t.85??7 

-.0002 

-.92E 

06 

.I79763B7E-0? 

1 .55360 

-.9003 

-.30E 

06 

-.179763871-0? 

! .55360 

-.0003 

-.10E 

06 

.9H9H601 OF-03 

.30S9M 

.0020 

.1«E 

06 

-.>(9 J86O1OE-03 

.39598 

.0020 

.lAE 

06 


0 QFfiSTTY .3095(IE-M 

Al«>HA 8t, 759970 PHin S. 50923? 


9«»05271or-l5 -. 

39O0OS6OE-JC, 

.10975191E-19 


TM4C 

pHtn 

ZFTA 

TC 


.20H73855E-02 

1.85189 

-.0003 

-.27F 

06 

-.20873PS5F-0? 

1 .85189 

-.0003 

-.?7E 

06 

.I7539366E-02 

1 .5560? 

-.0009 

-.?9E 

06 

-. I7539366E-0? 

1 .SShO? 

-.0009 

-.29F 

06 

.96087Q53F-01 

.90887 

.0027 

.I3E 

06 

-.9^0«7053E-03 

.90087 

.0027 

.13E 

06 


0 DFMSJTY .36100F-I1 

Al PHA fl?.0?07h0 PHtH S, 723057 


69705I86F-1S -. 

10960383E-19 

.1?301 

50UF-19 


ll^AR 

FPFf) 

ZFTA 

TC 


.?0'>1?331E-0? 

1 .851 9p 

-.0009 

-.20E 

06 

-.?00l?l3iE-n? 

1.85196 

-.0009 

-.?0E 

06 

.176(i9ni9F-C2 

1.55856 

-.0005 

-.20F 

06 

-. 1 7h090loF-02 

1.55856 

-.0005 

-.20E 

06 

,07q79o/4pp-.f,^ 

.92979 

.0033 

.lOE 

06 

-.97979O92E-07 

.92970 

.0053 

.lOE 

06 


58 


ALTITUDE 215,0 DENSITY ,«2150E-U 

^HT 150.J1S658 ALPHA 82,202089 PHIB 8,02T89T 

,520ai704E-17 ,97ia<1515F-|6 ,12195I06E-1« ,1223S6«8C-14 


REAL 

•;970»0251C*06 

-,97080251E-O6 

-,95960a76E-06 

-,95960«76E-06 

,19SO«073E-05 

,1930«073E-05 


IMAG 

,2095l88lE-02 

-,20951881E-02 

,l7869946E-02 

•,17669946E-0? 

,50179953E-03 

-,50179953E-OS 


FRED 2ETA 7C 

1 .85117— •; 0085 — ■;Ti rgi » 
1.85117 -.0005 -.16E 06 

1.56120 -.0005 '*-.17f 06 

1.56120 -.0005 -.17E 06 

.44336 ,0038 * • ,82E 05 

.44336 .0038 .82E 05 


ALTITUDE 210.0 DENSTTY ,«9300E-11 
PHI 148,378/478 ALPHA 82,566771 PHIP 6,302888 


95409791E-17 - 

.5R532293E-15 -,i 

66439909E-15 

.U884353E-14 


REAL 

imag 

FPEO 

■ ZETA 


— 

.11245645E-05 

.209927596-02 

1.85098 

-.0005 

-.14E 

06 

.1 1245645F-05 

-.20O92759E-02 

1,85098 

-.0005 

-,1«E 

06 

,115129966-05 

.17737296F-02 

1.56394 

-.0006 

-,14E 

06 

.11512996E-05 

-.17737296F-02 

1.56394 

-.0006 

-,14E 

06 

.22758641E-05 

.5272b675E-03 

,46490 

.0043 

,70E 

05 

.22758641F-05 

-.52726675E-03 

,46490 

.0043 

,70C 

or 


ALTITUDF 205,0 DENSITY ,578S0r-ll 
PHI 1/16,581409 ALPHA 8?,B<44995 PHIB 6,552577 
.OOOOOOOOE 00 .1 l79fel20E-1/4 -, 1 49880UE-1 /4 , 1 9073251 €-1 4 


REAL 

I NAG 

FREQ 

ZETA 

TC 


-.12418217E-05 

.21035236E-02 

1,85092 

-,0006 

-,13E 

06 

-.12418217E-05 

-.21035236E-0? 

1.85092 

-,0006 

-,13E 

06 

-.13774609E-05 

.17R06583E-02 

1,56682 

-.0008 

-.12E 

06 

-.137746096-05 

-.178O6503E-O2 

1.5668? 

-.0008 

-.12E 

06 

.26192826E-05 

.55676695E-03 

,48991 

,0047 

,6ir tT5 

.26192H26E-05 

-.55676693E-03 

,40991 

.0047 

,61E 

05 

ALTITUDE 200, 

0 DENSITY .680506 

-11 





PHI 1/4/4, 9473B/4 ALPHA 83,1 2/41 /43 PHIB 6,774560 
-,16479873E-16 -, 1 1 51 8564E- 1 4 -, 8604228aE-l 5 , 1 4378359E-1 4 


PEAL 

-.13166784F-05 

-,l3166784E-05 

-.16446367F-05 

-,16446367E-05 

,296l3l51E-05 

,?9613151F-0S 


I MAG 

,?1079621E-02 

-,210796216-02 

,17878071F-02 

-,17B78071E-02 

.59051743E-05 

-,59051745E-03 


FREQ 

1 .85101 

1.85101 
1 .56988 
1 .56988 

.51654 

.51854 


ZETA 

-.0006 

-.0006 

-.0009 

-.0009 

.0050 

.0050 


TC 

-,12E 06 
-,12E 06 
-,97E 05 
•,W 05 
.54E 05 
,54E 05 
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ALTITUDE t95,0 HFUiSTTY ,«OSOOE*|l 
9HI l<n,«72856 ALPHA PS. 905599 PHIB 6.9706IB 


•.S295979fcE-16 - 

,629S0095E-15 , 

13322876E-19 

.197179I9E-I8 


PEAL 

IMAG 

F9F0 

ZETA 

TC 


•,lS989783E-05 

.21126252E-0? 

I.85129 

-.0006 

:in 

06 

-.1S969783E-0B 

-.21128252E-02 

1.85129 

-.0006 

-.12E 

06 

-.196S1992E-0S 

.17952558E-02 

1.57S18 

-.0011 

-.81E 

05 

-.19651992E-05 

-.17952558E-02 

1.57318 

-.0011 

-.81E 

05 

.SS101275E-05 

,6290399OE-OS 

.55125 

.0055 

,a8F 

05 

,3S101?75E-('S 

-,8?903990E-03 

,55123 

.0055 

.98E 

05 


ALTITUnC 190,0 RfcN'SlTV ,q«750F-tl 
PHI l92,t79P16 ALPHA HS,6855«7 PHIB 7.1S9171 


.52091709E-17 -, 

1992R903E-15 -.( 

S6613381F-15 

.69590B91F-15 


REAL 

THAG 

FRFO 

ZFTA 

TC 


-, 1 3327127F-05 

,211 75901E-02 

1 .85177 

-.0006 

-,12E 

06 

-.t33?7l27F-f'5 

-.21 175901E-0? 

1 .85177 

-.0006 

-,12E 

06 

-,23365l50E-f5 

. |H''301 19E-0? 

1 ,5767? 

-.0013 

-,68E 

05 

-.2336S150F-05 

-.180301 19E-0? 

1.57672 

-.0013 

-.68E 

05 

,56692?77E-n5 

.87192999F-03 

.58760 

.0055 

,93E 

05 

.36692277E-0S 

-.67192999E-03 

.58769 

.0055 

,93E 

05 


AITITU'>E 185,0 DENSITY •llr'HSE-iO 
PHI iO(),09?tii9 ALPHA B'A.9«S350 PHIB 7,?91270 
-.50SON98U-1H ,971 «951SF-Ib -,230371 ?8E-1 0 , 25008657F- 1 9 


HEAL 

-,127?in«E-0B 
-,1?7?1 1 19F-05 
-,?H019529E-C5 
-,?«01«529E-l5 
,!J0755693E-''S 
,90735803E-/'S 


TvAr 

,212279708-02 
-,?l??797or-02 
,191 19317F-U2 
-,1Bn'»3l7F-9? 

,7????9aor-o3 

-,722?2990r-05 


FHFH 

ZETA 

TC 


1,85259 

-.0006 

-,I3F 

06 

1.85259 

-.0006 

-.13F 

06 

1,5808? 

-.0015 

-,57E 

05 

1 ,58082 

-.0015 

-.57E 

05 

.63028 

.0056 

.39E 

05 

.63028 

.0056 

.39E 

05 


Al TITUDF 1R0,0 nF^Sny ,133R5E-10 
PHI 139,992383 ALPHA R9,2Ro0ft1 PHIB 7,917952 


31699970E-17 -. 

16653395F-15 -.12351 

231E-19 

.12963099E-19 


PFAL 

T^*AG 

FPFO 

ZETA 

TC 


.117169328-05 

.?1?83895E-0? 

1 .85358 

-.0006 

-.1«E 

06 

, 11716932E-05 

-.21283895F-0? 

1.85358 

-.0006 

-,19E 

06 

.333729?Rp-(i5 

.18202819E-02 

1 .58526 

-.0018 

-.98E 

05 

.33372928F-05 

-,18?0?819E-0? 

1.58526 

-.0018 

-.98E 

05 

.95088H69E-05 

, 77f(6?786F-03 

.67635 

,0058 

.35E 

05 

,9508RR60E-05 

-.776627R6F-03 

.67635 

.0058 

.S5F 

05 


60 


• LTTTlIDf’ 175.0 DENSITY .1600OE-I0 
PHI IJD, 090055 ALPHA SO. 602675 PHIB 7.529699 


-.3069ao70E-17 

.15aO<J30«E-10 .12090009E-15 

,lS050936F-tO 


PEAL 

IMAG 

FREQ 

ZETA 

TC 


-.102705S6E-05 

.21300P53E-0? 

1.S5505 

-.0005 

-.15F. 

06 

-. 102705S6E-05 

-.2130aS53E-02 

1.S5505 

-.0005 

-.15E 

06 

-.39992176E-05 

,1S300785E-02 

1.59009 

-.0022 

-.OOF 

05 

-.39992176E-05 

-.1S3007B3E-O2 

1.59009 

-.0022 

-.OOE 

05 

.50262762F-05 

.S3901761E-O3 

.7295? 

.0060 

.32E 

05 

.50262762E-05 

-.B3O01761E-03 

.72952 

.0060 

.32E 

05 

ALTITUDE 170 

.0 DFURITY .1O250E 

-10 





PHI lJfl.3O0'’J0 ALOHA SO. 957136 PHIH 7.62S477 


-.S1532003E-16 

.3l9lPOi?F.i5 -.70900050E-15 

,eiS610S9E-15 


HEAI 

IvAn 

FRFQ 

ZETA 

TC 


-.aoi266l3F-)6 

.21012331F-O2 

1 .55706 

-.0000 

-.19F 

06 

-.B'll26fll3E-ii6 

-.21012331F-O2 

1.65706 

-.0000 

-.19F 

06 

-.081 3150BE-''5 

, 1 O'll 030CF-0? 

1.59669 

-.0026 

-.33E 

05 

-.'I813150HF-(’5 

-, I 1 03.inF-02 

1 ,59669 

-.0026 

-.33F 

05 

.56500?^9E-m5 

.91 12B906F-03 

.79035 

.0062 

.2SE 

05 

,56500??9f-''5 

-.9H2B906E-03 

.79035 

.0062 

.?SF. 

05 

ALTITLDE 165 

.0 DFNSITY .2515PE 

-to 





PHT 137.62^370 ALPHA 65.33P633 PHIB 7.7J272B 


-.M275703F-16 

.20261570H-1 0 -,?/J98(»018E-l 0 

.32169347E-10 


PFAL 

I mac 

FPFQ 

ZFTA 

TC 


-.6<10O515E-06 

.21 087309E-02 

1 .85969 

-.0003 

-.25E 

06 

-.63t095l5E-(.6 

-.21 '467309E-0? 

1.65969 

-.0003 

-.25F 

06 

-.57707166F-O5 

.1H531003F-02 

1.60383 

-.0031 

-.28F 

05 

-.S7747166F-''*: 

-.l«*'3K'03F-O2 

1 .60383 

-.0031 

-.28E 

05 

.600581 17E-05 

,9«>'139505F-03 

.85717 

.0065 

.25E 

05 

.600581 17F-05 

-.'’90395051-03 

.8571 7 

.0065 

.25F 

05 

ALTITUDE 160 

.0 DFriSITY .27950F 

-10 





PHI 137.007'4'iH ALPHA 8S.76f>019 PHIB 7.7S7501 


059701721-16 - 

.652256035-15 -, 

43021 142F-15 

.78270879E-15 


PEAL 

P'AG 

FPEQ 

ZETA 

TC 


.01 327615F-.I6 

.21573 )97F-i)2 

1.86320 

-.0002 

-.39E 

06 

.'41 3276’5E-16 

-.21573097F-0? 

1 .86320 

-.0002 

-.39E 

06 

.6930103OF-(m 

. 1 8667?q5E-0? 

1 .61227 

-.0037 

-.?3E 

05 

.69301030F-05 

-. I8667295E-0? 

1 .61227 

-.0037 

-.23F 

05 

.730341951-05 

. 1 0793185E-02 

.93219 

.0068 

.22E 

05 

. 73030 195F-05 

-. 1 079318BF-0? 

.93219 

.0068 

.22E 

05 


61 


ALTITUDF 155,0 DfNSITY .3aC50E-10 
PHI 156.5316«6 alpha AA.2653H7 PHIH 7,PS71<»a 


«B572?57F-16 - 

,Hl0789PaE-t5 -.929H1 

178E-15 

.12398R93E-14 


REAL 

IMAC 

FRED 

ZFTA 

TC 


,23O9O6«0E-Ofe 

.21675716E-0? 

1.86821 

-.0001 

-,69C 

06 

,250906BHF-06 

-.P1675716E-02 

1.86821 

-.0001 

-,69E 

06 

.B3002n7F-05 

.1BA25628E-02 

1.62256 

-.0044 

-.19E 

05 

.8340211 7E-05 

-.18H2562PE-0? 

1 ,62256 

-.0044 

-.19E 

05 

, 857 ine 5 F -05 

.1 IH18305F-02 

1,01861 

.0073 

.19E 

05 

.3S71HA5E-05 

-.1 181P305F-02 

1.01861 

.0073 

,19E 

05 

ALTITUDE 150 

,U DFHSITV .41600F 

-10 





PHI 136,01*05^^ ALPHA ftfe,f*35a«0 PHIR 7,9?a5«« 


-.17347235E-17 

,6?i<5(tfi45F-15 

10408341F-14 

.1213811 

9F-14 


RFAL 

TmaG 

FREQ 

ZETA 

TC 


-.2003612PE- )6 

.21H00972F-02 

1 .87509 

-.0001 

-.79E 

06 

-.20036! ?9F-r,6 

-.21 AOO972E-02 

1 .87509 

-.0001 

-.79E 

06 

-.996150|OE-«:5 

, 190061 79E-0? 

1 .63471 

-.0052 

-.16E 

05 

-.PP6150! 0F-n5 

-. 190061 79F-02 

1.6X971 

-.0052 

-. 1 6E 

05 

. 10161H62F-r.(j 

. . 12961691E-0? 

1 .1 10R3 

,0078 

.16E 

05 

. t0161862E-0/4 

-.12961691F-02 

1.11483 

,0078 

.16E 

05 

ALTITllDfc 105 

,0 OFNPITY .51200E-10 





PHI 135.701305 ALPHA «7. 505000 PHIB 7.9fl0At6 
.1hR26«1HF-15 ,1505q<(5ht..,q -.U102230F-15 . 1 60RR?7SE-1 0 


REAl 

-.S7007?<>3F-('ft 
-,570rt72biE-Oft 
1 1 73Ba?HF-'.0 
-.1 i73R0?RF-00 
. 1?30«500F-00 
. 1230P500t--n/^ 


U'AG 

.21963P6OF-0? 

-.?19f»3?<i9F-02 

,102l5Mq7F-0? 

-.17P15897F-0? 

.102#>oq3PE-0? 

102609TPF-0? 


FMED 

7FTA 

TC 


1.8851 1 

-.0003 

-.28E 

06 

1.8851 1 

-.0003 

-,28E 

06 

1 .64930 

-.0061 

-.14E 

05 

1.64931) 

-.0061 

-,14F 

05 

1 .2240? 

.0086 

.13E 

05 

1.22402 

,0086 

,13E 

05 


ALTTTline lOP.O HFHSITY .63000F-tO 
PHI 135,3hl6'l? ALPHA fle.?RRPP7 PHIR P, 060363 
,?5p«0??7E-l5 .057O6700F-15 707767 1 PF- 1 5 . P799R5R3E-15 


PEAl 

IPOOOaROF-nS 
190O0?6OF-ii5 
13?R->579F-o/J 
1 32n5579F-P'4 
. t51R9606F-oa 
.IRlflO^OHF-'O 


T'-’Af, 

,?2l«5ooiF-0? 
-.221R5001F-0? 
. 19057O2at-0? 

1 9057O20F-0? 
.1571H2I7E-0? 
-.1571PP17F-0? 


FREG 

ZFTA 

TC 


1 .9(1016 

-.0009 

-.84E 

05 

1 .90(^16 

-.0009 

-.84E 

05 

1 .66658 

-.0068 

-.12E 

05 

1.66658 

-.0068 

-.12E 

05 

1.34628 

.0097 

.lOE 

05 

1 .34628 

.0097 

.lOF 

05 


G2 


ALTITMOe t35,0 DENSITY ,7<»I50F*10 
PHI 155,0601 n ALP»M flO.Pljn?? PHTP S. 140166 
.216P404SE-16 -.156S1900E-14 ,PBSl 784?F*15 . 1 eo?3717C-l 4 

... # 


PEAL 

INAC 

EPED 

2ETA 

TC 


-,5464t688E-05 

,22511249F-02 

1.92424 

-*0024 

-,29E 

05 

-,54641684E-05 

-,225n2«9E-02 

1,92424 

-,0024 

-.29E 

05 

-,130921l4E-04 

,19741637F-02 

1,64735 

•,0066 

-,12E 

05 

-,lS092tl4E-<>4 

-,19741637E-02 

1,64735 

-,0066 

-*12E 

05 

,lS5562e3E-04 

,17340669E-02 

1,44213 

,0107 

• 46E 

04 

.1S556243E-04 

-,173406696-02 

1,46215 

,0107 

,46E 

04 


ALTITUDE 130,0 DENSITY ,99450E-10 
PHI 134,792731 ALPHA 90,286037 PHIP «, 244909 
,45536491E-!5 ,?t649349F-|4 -,3fi«*;7«06F-l 5 ,22461 730E-14 


PEAL 

IMAG 

FRED 

7ETA 

TC 


-,13454942E-rt4 

,230445946-02 

1.96553 

-.0054 

-.12E 

05 

-,13454942E-04 

-.2304459ME-0? 

1,96555 

-.0054 

-.12E 

05 

-,»0B29459E-f)6 

,?0M'6?77E-02 

1,71491 

-.0004 

-,20E 

06 

-,6n429459F-06 

-,?01('6277E-02 

1.71491 

-.0004 

•*20E 

06 

,142672376-04 

,190306056-02 

1.62317 

.0075 

.tie 

05 

,142672377-1)4 

-.190306056-0? 

1,62317 

.0075 

• HE 

05 

Al.TlTUnE 125 
PHI 134,551464 
,571230426-15 - 

,0 density ,1260(16-09 

ALPHA 91.526200 PHTfl 4.! 

,527355946-15 -,444049216-15 

J46674 

,745027356-15 


peal 

I mat. 

FPEO 

ZFTA 

TC 


-,2560092«E-04 

.240172«iF.-02 

2,04419 

-.0107 

-,6?E 

04 

-,256009246-04 

-,240172416-02 

2.04419 

-.0107 

-,62E 

04 

,447111146-04 

.212570096-02 

1.40756 

.0210 

.36E 

04 

,4471 1 1146-04 

-,21237onoE-0? 

1 ,40756 

.0210 

.366 

04 

-,191101946-04 

,100504936-0? 

1.69809 

-.0096 

-.8SE 

04 

-,191 I(M94F-'l4 

-, 19950493P-0? 

1 ,69409 

-.0096 

-.83F 

04 

ALTITUDE 120 
PHI 134,524505 
,5354691 76-15 

,0 DF'J.SITY ,162006-09 

ALPHA 92,972074 PHIR 8,579343 

.02941 1746-15 -. 444089? 1 E- 1 5 ,118513266-14 


PEAL 

INAG 

ERE (3 

ZFTA 

TC 


-.283l5927E-('4 

.25424062E-02 

2.19370 

-.0110 

-.56E 

04 

-.243159276-04 

-,254240626-0? 

2.19370 

-.0110 

-.56E 

04 

,504335406-04 

.227351136-02 

1.93100 

.0222 

.32E 

04 

,504335406-04 

-.227351136-0? 

1.93100 

.022? 

.32E 

04 

",?21 17653E-04 

.202H2672E-02 

1 ,72270 

-,010« 

-.726 

04 

-.221176536-04 

-.20242672E-02 

1.72270 

-.4109 

-.726 

04 


AITITUOE !|5,0 DENSITY ,210001-09 
PHI 11<1, 127320 alpha 9/J, 57313? PHIH S,S056«7 
,3209238«F-15 , I 502067PE-1« . 1 55«31 22E-1 « ,22009anC-ia 


REAL 

I MAG 

FREO 

ZETA 

TC 


20060273E-O9 

*?851998tE-02 

2, 0172/4 

-,0070 

-,79E 

09 

20060273E-09 

-.28519983F-02 

2,'4172« 

-.0070 

•,79E 

09 

00«95l2«E-09 

,2<M20t69E-02 

2.09033 

.0168 

,3RE 

09 

rtoa9512«E-oa 

-,?01?0189E-02 

2.09933 

.0168 

• 39E 

09 

?0<l30851E-0« 

.20995099E-02 

1,73708 

-.0100 

-,78E 

09 

2003a851E-O9 

-,20495099E-02 

1,73708 

-.0100 

-,78E 

0/1 

ALTITUDE 110 

,0 nfNSTTY .27800E 

• 1)9 





PHI lS3,P««n3 ALPHA 9fe.3?772‘» PHIB 9,222803 


-.10198132E-15 

.12628787E-19 , 

09000000E OO 

.12669a95E-19 


RFAL 

I^AG 

FRF9 

ZETA 

TC 


-,l?352992E-0« 

,31QO8*;?flP-03 

2.70835 

-.0039 

-.13F 

05 

-,1235?092F-oa 

-.3199852BF-0? 

2.70635 

-.0039 

-.13E 

05 

,303272?nF-0/j 

.25877953F-0? 

2.17173 

.0118 

.52C 

09 

.3O327220E-O/4 

-.25877/IR3F-02 

2.17173 

.0118 

,52E 

09 

-.179751786-09 

, .20891223E-02 

1.79578 

-.0087 

-,89E 

09 

-,17975178F-0a 

-.20891223E-O? 

1 .79578 

-,0087 

-.89F 

09 


ALTITUOP 105,0 OFNSITY ,368S0r-n9 
PHI 13S,7802«1 ALOMA ow.lAOSOH PHI« P.TSPOOP 
,l'>3'J2l87e-1« ,?3'l53aMF-i/| 22200460F -1 , 378^80 1 3E-1 0 


PEAL 

-,688905aHF-05 

-.OflSOOSaSF-''*? 

,22«55897F-00 

,2?/l55897E-'l/4 

-,15*;8880?E-na 

t5S888«2E-n/i 


THAG 

.3827UM89E-02 

-,i827««89F-02 

,278<l8P52F-0? 

-,?7898fl52F-n? 

.2O736389E-02 

-,2n736589e-o? 



FRFQ 

ZFTA 

TC 


3. 

,08157 

-,0019 

-.23E 

05 

3, 

,06157 

-.0019 

-.23E 

05 

2. 

,33337 

.0081 

• 71E 

09 

2. 

,33337 

,0081 

.Tlf 

0« 

1 . 

,75013 

-.0075 

-,10E 

05 

1. 

,75013 

-.0075 

-,10E 

05 


AUTITUDfc 100.0 re-NSTTY .SOlSOp-OP 
PHT t33.8377!l2 ALPHA 1 00,098579 PHIP 10,530307 
. l0g3703?F-l/| ,360R??0HF-15 -,???00O80E-15 , 1 1 729325E-1 0 


RFAL 

I MAR 

FHFQ 

ZETA 

TC 


285530 38F-r 5 

.9i538l5iF-n2 

3.99837 

-.0007 

-.58F 

05 

?ftS53038f.-05 

-.'M838151F-0? 

3.99837 

-.0007 

-,56F 

05 

18355717E-U/I 

.30237659E-02 

2.59669 

.0059 

.97F 

09 

16355717t-8a 

-.30237859E-02 

2.5988'I 

.0 059 

.97E 

09 

1350091 5E-09 

.20787831F.-02 

1.75077 

-.0065 

-.12E 

05 

1 350091 3F-09 

-.?,)787H3lF-n? 

1.75077 

-.0065 

-.12E 

05 


64 


ALTITUDE «*»,0 Op^SUV 

PHI 155.516696 ALPHA 101,97a562 PHIfl 11,659?59 


.59577716E-10 

.22209/160F-19 -,66617*«2E-15 

.96291970E-19 


REAL 

ihag 

FREQ 

ZETA 

TC 


.3666a515E-06 

.961«869PE-02 

9.ft«65'l 

.0001 

.95E 

06 

,56B««395E-06 

-,9R196699E-02 

9.09659 

.0001 

.95E 

06 

.115249««E-0/l 

.5571379iJE-02 

2.65559 

.0039 

.19E 

05 

.1152«96ttE-ifa 

-.33715799E-02 

2.65.359 

.0039 

.19E 

05 

-.1189592BE-09 

.2079559RF-02 

1.79753 

-.0057 

-.15E 

05 

-,llB93P26E-09 

-.20793396E-02 

1.79755 

-.0057 

-.15E 

05 

ALTITUDE 90 

.0 DENSITY ,10095f-06 




PHI 133.920090 

ALPHA lOl.TOBRBO 

PHIH 13.371772 



.6552917PE-19 

.566213790-19 -. 2??09960E-1 5 

.66651509E-19 


REAL 

IKAG 

FREQ 

ZFTA 

TC 


.3? 19951 OF -OR 

.5671O3O5F-0? 

9.75590 

.0006 

.50E 

05 

,32l«951BE-oR 

-.567103050-02 

9.75590 

.0006 

.50E 

05 

.76309357E-'i5 

. 4B'I71P90E-0? 

3.22693 

.0020 

.21F 

05 

.76509357F-H5 

-.3B97)99nF-02 

3.22693 

.0020 

.21F 

05 

10R«9R«HF-ro 

,2'»7‘3n??E-0? 

1.75061 

-.9052 

-.15E 

05 

-.lOBaOBAftE-DU 

-.20751 122F-0? 

1.73661 

-.0052 

-.15F 

05 

ALTITUDE RS 

.0 OFNSITY .150500-00 





PHI 153,35?6?0 ALPHA l(’5.3'i2«3« PHIR 16,1/10770 
.173M7355E-1R .P3R69795E-1 « .2PP0OO60E-1 R ,?«03553lE-ia 


real 

,61567?39F-05 

.6l567230£-'5 

,a3«06«H9E-rs 

,'UB066R9t-<i5 

-,t 0*5370l3E-'''» 
-.10537«l3L-('U 


l*'*Ar, 

.hBPoninF.-oa 

-./>6?9138or-02 

.051RP01/IE-0? 

“.051HPP1/1E-0? 

.?052?P70fc-02 

•,209??P7«F-oa 



FRf Q 

7FTA 

TC 


I. 

.71507 

,0009 

.26E 

05 

5, 

,71507 

.0009 

.26E 

05 

3. 

,70179 

.0010 

.36E 

05 

3. 

.70179 

.0910 

.36E 

05 

1, 

,71759 

-.0051 

-.15E 

05 

1. 

,71759 

-.0051 

-.15E 

05 


OrtiSTTY ,2«00«)E-OB 
ALPHA lOh. 76/1001 PHIR 21.176511 
-.20P61S70F-1U -.17763R6BE-10 ,56312262E-1« 


ALTITline BO.O 
PHI 133.3003PO 
.P723515PE-1/1 


r-'Ao 

.H0Pf'5Rl3F-n? 

-.RaPRRSlSF-O? 

.551RH136F-0? 
-,55l0H136fc-O2 
, tP«66PS6F-P? 
-,19B66RS6t-02 


FREQ 

ZFTA 

TC 


7.99799 

.001? 

.16E 

05 

7.09709 

.0012 

.16E 

05 

9.60069 

.0095 

.ME 

06 

9.60060 

.0093 

.11F 

06 

1.65906 

-.0050 

-.1«E 

05 

1 .65900 

-.0056 

-,19F 

05 


PEAL 

. 1005022RF-'<I 
.1005iJ22RF-;a 
.13R77365L-U5 
• 13B77365E-''*; 

-.1 19ai965F-;)0 
-,1 1««lO65fc-i.0 


65 













APPROVAL 


EQUILIBRIUM AND STABILITY OF A SATELLITE INFLUENCED 
BY GRAVITATIONAL AND AERODYNAMIC TORQUES 

By Zachary J. Galaboff 


The information in this I'eport has been reviewed for technical 
content. Review of any information concerning: Department of Defense 
or nuclear energy activities or programs has been made by the MSFC 
Security Classification Officer. This report, in its entirety, has been 
determined to be unclassified. 



Director, Systems Dynamics Labox’atory 


74 


■<>^U.S. GOVERNMENT PRINTING OFFICE 1981—740-066/316 REGION N0.4 


